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The modified newtonian dynamics (MOND) paradigm of Milgrom can boast of a number of 
successful predictions regarding galactic dynamics; these are made without the assumption that 
dark matter plays a significant role. MOND requires gravitation to depart from Newtonian theory 
in the extragalactic regime where dynamical accelerations are small. So far relativistic gravitation 
theories proposed to underpin MOND have either clashed with the post-Newtonian tests of general 
relativity, or failed to provide significant gravitational lensing, or violated hallowed principles by 
exhibiting superluminal scalar waves or an a priori vector field. We develop a relativistic MOND 
l/^ , inspired theory which resolves these problems. In it gravitation is mediated by metric, a scalar 

■ field and a 4-vector field, all three dynamical. For a simple choice of its free function, the theory 
' has a Newtonian limit for nonrelativistic dynamics with significant acceleration, but a MOND limit 
I when accelerations are small. We calculate the /3 and 7 PPN coefficients showing them to agree 

^jQ. with solar system measurements. The gravitational light defiection by nonrelativistic systems is 

r—{ ' governed by the same potential responsible for dynamics of particles. To the extent that MOND 

I successfully describes dynamics of a system, the new theory's predictions for lensing by that system's 

, visible matter will agree as well with observations as general relativity's predictions made with a 

' dynamically successful dark halo model. Cosmological models based on the theory are quite similar 

^N) ' to those based on general relativity; they predict slow evolution of the scalar field. For a range 

, of initial conditions, this last result makes it easy to rule out superluminal propagation of metric, 

■ scalar and vector waves. 
> ■ 

' PACS numbers: 95.35.+d,95.30.Sf, 98.62.Sb, 04.80.Cc 
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Tj- , I. INTRODUCTION 

O ; 

. In the extragalactic regime, where Newtonian gravitational theory would have been expected to be an excellent 
Qh' description, accelerations of stars and gas, as estimated from Doppler velocities and geometric considerations, are as a 
Q , rule much larger than those due to the Newtonian field generated by the visible matter in the system This is the 

■ "missing mass" problem or "acceleration discrepancy" 0. It is fashionable to infer from it the existence of much 
[ dark matter in systems ranging from dwarf spheroidal galaxies with masses IO^Mq to great clusters of galaxies in 

d the regime HQ. And again, galaxies and clusters of galaxies are found to gravitationally lense background 

^ sources. When interpreted within general relativity (GR), this lensing is anomalously large unless one assumes the 
presence of dark matter in quantities and with distribution similar to those required to explain the accelerations of 
stars and gas. Thus extragalactic lensing has naturally been regarded as confirming the presence of the dark matter 
suggested by the dynamics. 

But the putative dark matter has never been identified despite much experimental and observational effort 0. 
This raises the possibility that the acceleration discrepancy as well as the gravitational lensing anomaly may reflect 
departures from Newtonian gravity and GR on galactic and larger scales. Now alternatives to GR are traditionally 
required to possess a Newtonian limit for small velocities and potentials; thus the acceleration discrepancy also raises 
the possibility that the correct relativistic gravitational theory may be of a kind not generally considered hitherto. 

In the last two decades Milgrom's modified Newtonian dynamics (MOND) paradigm has gained recognition 

as a successful scheme for unifying much of extragalactic dynamics phenomenology without invoking"dark matter" . 
In contrast with earlier suggested modifications of Newton's law of universal gravitation [HIIll^JlJ, MOND 
is characterized by an acceleration scale ao, not a distance scale, and its departure from Newtonian predictions is 
acceleration dependent: 

/i(|a|/ao)a = -V$N- (1) 
Here $n is the usual Newtonian potential of the visible matter, while fl{x) « x for x ^ 1 and fl{x) — )■ 1 for x ^ 1. 
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Milgrom estimated ao « 1 x 10~* cm s~^ from the empirical data. In the laboratory and the solar system where 
accelerations are strong compared to Oq, formula goes over to the Newtonian law a — — V$n- 

Milgrom constructed formula to agree with the fact that rotation curves of disk galaxies become flat outside their 
central parts. That far out a galaxy of mass M exhibits an approximately spherical Newtonian potential. The scales 
are such that |V$n| ~ GMr^^ <C ao in this region, and so Eq. ^ with fi{x) « x gives |a| « (GAf ao)^^^r~^ which has 
the r dependence appropriate for the centripetal acceleration /r of a radius independent rotational velocity Vc — an 
asymptotically flat rotation curve. In fact one obtains the relation M = (Gao)~^Wc^ which leads to the prediction 
that for any class of galaxies with a constant mass to luminosity ratio T in a specified spectral band, the lumino sity 
in that band should scale as Vc'^. And indeed, there exists an empirical law of just this form: the Tully-Fisher law 
(TFL) relating near infrared (H-band) luminosity Ljj of a spiral disk galaxy to its rotation velocity, Ljj oc Vc'^, with 
the proportionality factor being constant within each galactic morphology class. 

This version of the TFL was established only after MOND was enunciated . It is in harmony with the MOND 
prediction in two ways. First, the infrared light of a galaxy comes mostly from cool dwarf stars which make up most of 
its mass (hence giving a tight correlation between M of the predicted relation and Lh of the empirical law). Second, 
the proportionality coefhcient varies from class to class as would be expected from the observed correlation between 
T of a galaxy and its morphology. 

In the alternative dark matter paradigm (which casts no doubt on standard gravity theory), flat rotation curves 
are explained by assuming that every disk galaxy is nested inside a roundish spherical halo of dark matter 16] whose 
mass density drops approximately like r~^. The halo is supposed to dominate the gravitational field in the outer 
parts of the galaxy. This makes the Newtonian potential approximately logarithmic with radius in those regions, 
thereby leading to an asymptotically flat rotation curve. In practice the dark halo resolution works only after fine 
tuning. It is an observational fact that for bright spiral galaxies the rotation curve in the optically bright region is 
well explained in Newtonian gravity by the observed matter |l7j| . But, as mentioned, in the outer regions the visible 
matter's contribution must be dwarfed by the halo's. So fine tuning is needed between the dark halo parameters 
(velocity dispersion and core radius) and the visible disk ones 0, 0| . 

This fine tuning problem is exacerbated by the TFL Lh oc Vc'^. Because the infrared luminosity comes from the 
visible matter in the galaxy, but the rotation velocity is mostly set by the halo, the TFL also requires fine tuning 
between halo and disk parameters. The standard dark matter explanation of the r^^ profile of an halo is that it arises 
naturally from primordial cosmological perturbations pol | . The visible galaxy is regarded as forming by dissipational 
collapse of gas into the potential well of the halo. The fine tunin g in entioned is then viewed as resulting from the 
adjustment of the halo to the gravitation of the incipient disk [l9l I2H. But the TFL is observationally a very sharp 
correlation; in fact, it is the basis for one of the most reliable methods for gauging distances to spiral galaxies. 
Such sharpness is hardly to be expected from statistical processes of the kind envisaged in galaxy formation, a point 
emphasized by Sanders |22| . So in the dark matter picture the TFL is something of a mystery. 

There are other MOND successes. Milgrom predicted early that in galaxies with surface mass density well below 
aoG^^, the acceleration discrepancy should be especially large In dwarf spirals this property was established em- 
pirically years later l23| , and it is now known to be exhibited by a large number of low surface brightness galaxies |24j . 
Another example: MOND successfully predicts the detailed shape of a rotation curve from the observed matter (stars 
and gas) distribution on the basis of a single free parameter, T, down to correlating features in the velocity field with 
those seen in the light distribution [25I Eg . . This is especially true in the case of low surface mass density 

disk galaxies for which MOND's predictions are independent of the specific choice of fi{x) ,£9,], and these MOND 
theoretical rotation curves fit the observed curves of a number of low surface brightness dwarf galaxies [l^, HO, |3l| 
very well. By contrast, the dark halo paradigm requires one or two free parameters apart from T to approximate 
the success of the MOND predictions |33|. In fact, even when the empirical data is analyzed within the dark halo 
paradigm, it displays the preferred acceleration scale ao of MOND . 

Occasionally doubt has been cast on MOND's ability to describe clusters of galaxies properly [s^]. Many of these 
exhibit accelerations not small on scale ao, yet conventional analysis suggests they contain much dark matter in 
opposition to what MOND would suggest. Sanders has recently reanalyzed the problem "s^ with the conclusion 
that these clusters may contain much as yet undiscovered baryonic matter in the core which should be classed as 
"visible" in connection with MOND. Other MOND successes, outside the province of disk galaxies, have been reviewed 
elsewhere HHHIs^. 

So the simple MOND formula ^ is very successful. But it is not a theory. Literally taken the MOND recipe 
for acceleration violates the conservation of momentum (and of energy and of angular momentum) ■ And MOND 
entails a paradox: why does the center of mass of a star orbit in its galaxy with anomalously large acceleration 
given by Eq. with /i <C 1, while each parcel of gas composing it is subject to such high acceleration that is 
should, by the same formula, be accelerated Newtonially ? In short, the MOND formula is not a consistent 
theoretical scheme. Neither is MOND, as initially stated, complete. For example, it does not specify how to calculate 
gravitational lensing by galaxies and clusters of galaxies. As is well known, in standard gravity theory light deflection 
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is well described only by relativistic theory (GR). And whereas Newtonian cosmological models work well for part of 
the cosmological evolution, MOND cosmological models built in analogy with their Newtonian counterparts, though 
sometimes agreeing with phenomenology Is.^ , can yield peculiar predictions 1^3 (but see Ref.ls^. In short, a complete, 
consistent theoretical underpinning of the MOND paradigm which accords with observed facts, and is also relativistic, 
has been lacking. 

This lack is being resolved in measured steps. A first step was the lagrangian reformulation of MOND [s^l called 
AQUAL (see Sec. Ill A\ . AQUAL cures the nonconservation problems and resolves the paradox of the galactic motion 
of an object whose parts accelerate strongly relative to one another; it does so in accordance with a conjecture of 
Milgrom 0- And for systems with high symmetry AQUAL reduces exactly to the MOND formula 

A relativistic generalization of AQUAL is easy to construct with help of a scalar field which together with the metric 
describes gravity [s^ (see Sec. Ill CTI below). It reduces to MOND approximately in the weak acceleration regime, to 
Newtonian gravity for strong accelerations, and can be made consistent with the post-Newtonian solar system tests 
for GR. But relativistic AQUAL is acausal: waves of the scalar field can propagate superluminally in the MOND 
regime (see the appendix of Ref. [^^ or AppendixElhere). The problem can be traced to the aquadratic kinetic part of 
the lagrangian of the theory which mimics that in the original AQUAL. A theory involving a second scalar field, PCG, 
was thus developed to by pas s the problem 40 400.^] (see Sec. Ill 0*21 below) . PCG may be better behaved causally 
than relativistic AQUAL [ij, but it brings woes of its own. It is marginally in conflict with the observed perihelion 
precession of Mercury and in common with relativistic AQUAL, PCG predicts extragalactic gravitational lensing 
which is too weak if there is indeed no dark matter. This last problem is traceable to a feature common to PCG and 
relativistic AQUAL: the physical metric is conformal to the metric appearing in the Einstein-Hilbert action [43| . 

One way to sidestep this problem without discarding the MOND features is to exploit the direction defined by the 
gradient of the first scalar field to relate the physical metric to the Einstein metric by a disformal transformation (see 
Ref. ji^l or Sec. Ill C"3l below). But it turns out that with this relation the requirement of causal propagation acts 
to depress gravitational lensing |44j . rather than enhancing it as is observationally required. The persistence of the 
lensing problem in modified gravitational theories has engendered a folk theorem to the effect that it is impossible for 
a relativistic theory to simultaneously incorporate the MOND dynamics, observed gravitational lensing and correct 
post-Newtonian behavior without calling on dark matter ,45„ ,46b 4L JS] ■ 

Needless to say, this theorem cannot be proved ■ Indeed, by the simple device of relating t he p hysical and Einstein 
metrics via a disformal transformation based on a constant time directed 4-vector, Sanders l5C| has constructed an 
AQUAL like "stratified" relativistic theory which gives the correct lensing while ostensibly retaining the MOND 
phenomenology and consistency with the post-Newtonian tests. Admittedly Sanders' stratified theory is a preferred 
frame theory, and as such outside the traditional framework for gravitational theories. But it does point out a trail 
to further progress. 

The present paper introduces TeVe S, a new relativistic gravitational theory devoid of a priori fields, whose nonrel- 
ativistic weak acceleration limit accords with MOND while its nonrelativistic strong acceleration regime is Newtonian. 
TeVe S is based on a metric, and dynamic scalar and 4-vector fields (one each) ; it naturally involves one free function, 
a length scale, and two positive dimensionless parameters, k and K. TeVe S passes the usual solar system tests of GR, 
predicts gravitational lensing in agreement with the observations (without requiring dark matter), does not exhibit 
superluminal propagation, and provides a specific formalism for constructing cosmological models. 

In Sec. ^ we summarize the foundations on which a workable relativistic formulation of MOND must stand. We 
follow this with a brief critical review of relativistic AQUAL, PCG and disformal metric theories, some of whose 
elements we borrow. Sec. IIII Al builds the action for TeVe S while Sec. IIII Bl derives the equations for the metric, 
scalar and vector fields. In Sec. IIII CI we demonstrate that TeVe S has a GR limit for a range of small k and K. 
This is shown explicitly for cosmology fSec. IIII C l|l and for quasistatic situations like galaxies fSec. IIII C 2|) . All the 
above applies for any choice of the free function; in Sec. IIII El we make a simple choice for it which facilitates further 
elaboration. For spherically symmetric systems the nonrelativistic MOND limit is derived in Sec. IIVBI while the 
Newtonian limit is recovered for modestly small k in Sec. II V CI The above conclusions are extended to nonspherical 
systems in Sec. II V Dl Sec. |V] shows that the theory passes the usual post-Newtonian solar system tests. Sec. IVII 
demonstrates that for given dynamics, TeVeS gives the same gravitational lensing as does a dynamically successful 
dark halo model within GR. In Sec. lVIII we discuss TeVe S cosmological models with flat spaces showing that they are 
very similar to the corresponding GR models (apart from the question of cosmological dark matter which is left open), 
and demonstrating that the scalar field evolves little, and so can be taken to be small and positive. As discussed next 
in Sec. IVIlTI this last conclusion serves to rule out superluminal propagation in TeVeS. 
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II. THEORETICAL FOUNDATIONS FOR THE MOND PARADIGM 



A. AQUAL: nonrelativistic field reformulation of MOND 



However successful empirically when describing motions of test particles e.g. stars in the collective field of a 
galaxy, formula ^ is not fully correct. It is easily checked that a pair of particles accelerating one in the field of 
the other according to |^ does not conserve momentum. Thus the MOND formula by itself is not a theory. It 
is, however, a simple matter to construct a fully satisfactory nonrelativistic theory for MOND ( |39j|). Suppose we 
retain the Galilean and rotational invariance of the Lagrangian density which gives Poisson's equation, but drop the 
requirement of linearity of the equation. Then we come up with 

Here p is the mass density, ao is a scale of acceleration introduced for dimensional consistency, and / is some function. 
Newtonian theory (Poisson's equation) corresponds to the choice f{y) = y. From Eq. (O follows the gravitational 
field equation 

V • [/i(|V$|/ao) V*] = 47rGp, (3) 

where il{\/y) = df{y)/dy. Because of its AQUAdratic Lagrangian, the theory has been called AQUAL 4]. The form 
of / and the value of ao must be supplied by phenomenology. We assume 

f{y)^\ lym (4) 



For systems with spherical, cylindrical or planar geometry, Eq. © can be integrated once immediately. With the 
usual prescription for the acceleration, 

a=-V4>, (5) 

the solution corresponds precisely to the MOND formula This is no longer true for lower symmetry. However, 
numerical integration reveals that is approximately true, in most cases to respectable accuracy |5lj . 

The mentioned inexactness of Eq. |^ for systems such as a discrete collection of particles is at the root of the 
mentioned violation of the conservation laws. Because AQUAL starts from a Lagrangian, it respects all the usual 
conservation laws (energy, momentum and angular momentum), as can be checked directly '39^. This supplies the 
appropriate perspective for the mentioned failings of MOND. AQUAL also supplies the tools for showing that New- 
tonian behavior of the constituents of a large body, e.g. a star, is consistent with non-Newtonian dynamics of the 
latter's center of mass in the weak collective field of a larger system, e.g. a galaxy. 

To summarize, whenever parts of a system devoid of high symmetry move with accelerations weak on scale Qq, the 
field which defines their accelerations is to be calculated by solving the AQUAL equation Q. AQUAL then 
becomes the nonrelativistic field theory on which to model the relativistic formulation of the MOND paradigm. 



B. Principles for relativistic MOND 



A relativistic MOND theory seems essential if gravitational lensing by extragalactic systems and cosmology are to 
be understood without reliance on dark matter. What principles should the relativistic embodiment of the MOND 
paradigm adhere to ? The following list is culled from those suggested by Bekenstein 0, 0|, Sanders and 
Romatka [53l |. 



1, Principles 



• Action principle The theory must be derivable from an action principle. This is the only way known to guarantee 
that the necessary conservation laws of energy, linear and angular momentum are incorporated automatically It is 
simplest to take the action as an integral over a local lagrangian density. A nonlocal action has been tried |47l |. but 
the resulting theory fails on account of gravitational lensing. 
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• Relativistic invariance Innumerable elementary particle experiments provide direct evidence for the universal 
validity of special relativity. The action should thus be a relativistic scalar so that all equations of the theory are 
relativistically invariant. Implied in this is the correspondence of the theory with special relativity when gravitation 
is negligible. This proviso rules out preferred frame theories. 

• Equivalence principle As demonstrated with great accuracy (1 part in 10^^) by the Eotvos-Dicke experi- 
ments free particles with negligible self-gravity fall in a gravitational field along universal trajectories (weak 
equivalence principle). For slow motion (the case tested by the experiments), the equation a = — V^, which encap- 
sulates the universality, is equivalent to the geodesic equation in a (curved) metric gap with gu ~ —1 — 2$. For light 
propagating in a static gravitational field, such a metric would predict that all frequencies as measured with respect to 
(w.r.t.) observers at rest in the field undergo a redshift measured by This is experimentally verified |55| to 1 part 
in 10^. It thus appears that a curved metric gap describes those properties of spacetime in the presence of gravitation 
that are sensed by material objects. According to Schiff's conjecture [s^. Is^. this implies that the theory must be 
a metric theory, i.e., that in order to account sfor the effects of gravitation, all nongravitational laws of physics, e.g. 
electromagnetism, weak interactions, etc. must be expressed in their usual laboratory forms but with the metric gap 
replacing the Lorentz metric. This is the Einstein equivalence principle '53|. 

• Causality So as not to violate causality and thereby compromise the logical consistency of the theory, the 
equations deriving from the action should not permit superluminal propagation of any measurable field or of energy 
and linear and angular momenta. Superluminal here means exceeding the speed which is invariant under the Lorentz 
transformations. By Lorentz invariance of Maxwell's equations this is also the speed of light. In curved space, where 
curvature can cause waves to develop tails, the maximal speed is that of wavefronts, typically that of the high frequency 
components. 

• Positivity of energy Fields in the theory should never carry negative energy. From the quantum point of view 
this is a precaution against instability of the vacuum. This principle is usually taken to mean that the energy density 
of each field should be nonnegative at each event (local positivity) . The fact that the gravitational field itself cannot 
be generically assigned an energy density shows that this popular conception is overly stringent. A more useful 
statement of positivity of energy is that any bounded system must have positive energy (global positivity instead of 
the stronger local positivity). For example, the gravitational field can carry negative energy density locally (at least 
in the Newtonian conception), yet for pure gravity and in some cases in the presence of matter, a complete gravitating 
system is subject to the positive energy theorems [s^]. Also, there are examples of scalar fields whose local energy 
density is of indefinite sign, yet a complete stationary system of such fields with sources has positive mass [s^ . Of 
course, local positivity implies global positivity. 

• Departures from Newtonian gravity The theory should exhibit a preferred scale of acceleration below which 
departures from Newtonian gravity should set in, even at low velocities. 

2. Requirements 

The relativistic embodiment of MOND should predict a number of well established phenomena. For example, we 
expect the following. 

* Agreement with the extragalactic phenomenology: The nonrelativistic limit of the theory should make predictions 
in agreement with those of the AQUAL equation, which is known to subsume much extragalactic phenomenology. 
This is checked for TeVe S in Sec. HVBl 

★ Agreement with phenomenology of gravitational lenses: The theory should predict correctly the lensing of electro- 
magnetic radiation by extragalactic structures which is responsible for gravitational lenses and arcs. In particular, it 
should give predictions similar to those of GR within the dark matter paradigm. This point is established for TeVe S 
in Sec. |Vl| 

* Concordance with the solar system: The theory should make predictions in agreement with the various solar 
system tests of relativity [H^: deflection of light rays, time delay of radar signals, precessions of the perihelia of the 
inner planets, the absence of the Nordtvedt effect in the lunar orbit, the nuUness of aether drift, etc. TeVeS is 
confronted with the first three tests in Sec. 

* Concordance with binary pulsar tests : The theory should make predictions in harmony with the observed pulse 
times of arrival from the various binary pulsars. These contain information about relativistic time delay, periastron 
precession and the orbit's decay due to gravitational radiation. They thus constitute a test of the strong potential 
limit of the theory. 

★ Harmony with cosmological facts: The theory should give a picture of cosmology in harmony with basic empirical 
facts such as the Hubble expansion, its timescales for various eras, existence of the microwave background, light 
element abundances from primordial nucleosynthesis, etc. The similarity of cosmological evolution in GR and in 
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TeVe S is established in Sec. IVIII though the problem of how to eliminate cosmological dark matter with TeVe S is 
left open. 

C. Some antecedent relativistic theories 

It is now in order to briefly review some of the previous attempts to give a relativistic theory of MOND. This will 
introduce the concepts to be borrowed by TeVe S, and help to establish the notation and conventions that we shall 
follow. A metric signature +2, and units with c = 1 are used throughout this paper. Greek indeces run over four 
coordinates while Latin ones run over the spatial coordinates alone. 

1. Relativistic AQUAL 

It is well known that theories constructed, for example, by using a local function of the scalar curvature as Lagrangian 
density, have a purely Newtonian limit for weak potentials. So if we steer away from nonlocal actions, then AQUAL 
behavior cannot arise from merely modifying the gravitational action. The theory one seeks has to involve degrees of 
freedom other than the metric. 

In the first relativistic theory with MOND aspirations, relativistic AQUAL the physical metric was taken 
as conformal to a primitive (Einstein) metric gap, i-e., gap — gap with il) a real scalar field. In order not to break 
violently with GR, which is well tested in the solar system (and to some extent in cosmology), the gravitational action 
was taken as the Einstein-Hilbert's one built out of gap- The MOND phenomenology was implanted by taking for 
the Lagrangian density for ■0 

'Cv,=-^^/(iVV',aV',/3), (6) 

where / is some function (not known a priori)^ and L is a constant with dimensions of length introduced for dimensional 
consistency. Note that when /(y) — y, is just the lagrangian density for a linear scalar field, but in general is 
aquadratic. 

To implement the universality of free fall, one must write all lagrangians of matter fields using a single metric, 
which is taken as gap (not gap which choice would make the theory GR). Thus for example, the action for a particle 
of mass m is taken as 

Sm = -m I {-gapdx'^dxPf'^ (7) 



Hence test particle motion is nongeodesic w.r.t. gap but, of course, geodesic w.r.t. gap- Evidently this last is the 
metric measured by clocks and rods, hence the physical metric. Addition of a constant to -0 merely multiplies all 
masses by a constant (irrelevant global redefinition of units), so that the theory is insensitive to the choice of zero of 
^- 

For slow motion in a quasistatic situation with nearly flat metric gap, and in a weak field "0, {— gapdx" dx^Y^^ « 
(1 + $N + "0 ~ v^/2)(it, were — —{gtt + l)/2 is the Newtonian potential determined by the mass density p through 
the linearized Einstein equations for gap, and v is the velocity defined w.r.t. the Minkowski metric which is close to 
gap- Thus the particle's lagrangian is m(v^/2 — $n — 0); this leads to the equation of motion 

aw-V($N + V^)- (8) 

How is ip determined ? For stationary weak fields the Lagrangian density for ip, including a point source of physical 
mass M at r = 0, is from the above discussion and Eqs. I©- HI), 

= {L\V^f) - ^MS{r). (9) 

Comparing Eqs. Q and l(2Jl we conclude that ip here corresponds to $ of mass M as computed from AQUAL's 
Eq. |j2Jl, provided we take f = f and L — l/oo. Whenever I'Vipl 3> |V$n| {^n is the Newtonian potential of the same 
mass distribution), the equation of motion ^ reduces to jSJ, and we obtain MOND like dynamics. For the choice of 
MOND function |0J the said strong inequality is automatic in the deep MOND regime, IVf/'l ^ O-o, because jl <^ 1 
there. 
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In the regime \'Vil'\ ^ ag, /i « 1 and f{y) « ?/ so that ■0 reduces to <i>N- It would seem from Eq. ^ that a particle's 
acceleration is then twice the correct Newtonian value. However, this just means that the measurable Newton's 
constant Gn is twice the bare G appearing in £^ or in Einstein's equations. It is thus clear, regarding dynamics, that 
the relativistic AQUAL theory has the appropriate MOND and Newtonian limits depending on the strength of 'Vtp. 

But relativistic AQUAL has problems. Early on ,4 .39». (4"3 | it was realized that ip waves can propagate faster than 
light. This acausal behavior can be traced to the aquadratic form of the lagrangian, as explained in Appendix fXI 
A second problem 113 issues from the conformal relation ga/3 = e^'^gap- Light propagates on the null cones of 
the physical metric; by the conformal relation these coincide with the lightcones of the Einstein metric. This last 
is calculated from Einstein's equations with the visible matter and -0 field as sources. Thus so long as the field 
contributes comparatively little to the energy momentum tensor, it cannot affect light deflection, which will thus be 
that due to the visible matter alone. But in reality galaxies and clusters of galaxies are observed to deflect light 
stronger than the visible mass in them would suggest. Thus relativistic AQUAL fails to accurately describe light 
deflection in situations in which GR requires dark matter. It is thus empirically falsified. 

Relativistic AQUAL bequeaths to TeVe S the use of a scalar field to connect Einstein and physical metrics, a field 
which satisfies an equation reminiscent of the nonrelativistic AQUAL Eq. iPJ . 



2. Phase coupling gravitation 

The Phase Coupled Gravity (PCG) theory was proposed 0, l40l li^ in order to resolve relativistic AQUAL's 
acausality problem. It retains the two metrics related by gap = e^"^ gap, but envisages ?A as one of a pair of mutually 
coupled real scalar fields with the Lagrangian density (our definitions here differ slightly from those in Ref. ^ 

C^.A = -\ [5"''(AaA/3 + ??~'A2^,„VA/3) + V(A2)] (10) 

Here 77 is a real parameter and V a real valued function. The coupling between A and "0 is designed to bring about 
AQUAL-like features for small The theory receives its name because matter is coupled to ■0, which is proportional 
to the phase of the self- interacting complex field x = Ae^^l'''^ . 

Variation oi C^^a w.r.t. A leads to (all covariant derivatives and index raising w.r.t. gap) 

A'";„-ry-2AV;,„0'"- AV'(A2) = (11) 

In the variation w.r.t. we must include the Lagrangian density of a source, say a point mass M at rest at r = 
[c.f. S„, in Eq. (0)]: 

(A2g"/5V,/3)^„ = ^'e'/'M5(r) (12) 

The connection with AQUAL is now clear. For sufficiently small \r]\ the A'";q, term in Eq. (|ll|l becomes negligible, 
and the other two establish an algebraic relation between V'.qV''" . Substituting this in Eq. H12|l gives the 

AQUAL type of equation for ^0 that would derive from in Eq. 

The PCG Lagrangian's advantage over that of the relativistic AQUAL's is precisely in that it involves first deriva- 
tives only in quadratic form. This would seem to rule out the superluminality generating A"" dependent terms 
discussed in Appendix ^ In practice things are more complicated. A detailed local analysis employing the eikonal 
approximation 42] shows that there are superluminal -0 perturbations, for example when V" < 0. However, the same 
analysis shows that such superluminality occurs only when the background solution is itself locally unstable. This 
makes the said causality violation moot. 

One way to obtain MOND phenomenology from PCG is to choose V(^^) — — ie~^A^ with e a constant with 
dimension of energy. Although with this choice V" < which makes for unstable backgrounds, we only need this 
form for small A\ V can take different form for large argument. Then in a static situation with nearly flat gap and 
weak -0, Eqs. lfTT)) - ((T^ reduce to 

V^^ - 7r^^(Vi/')^ + e^^A^ = 0, (13) 
V ■{A^V'il))^T]^M&{v). (14) 

The spherically symmetric solution of Eqs. H13() - (|14|l is 

A = (Ke/r)i/2; d^A/dr = {^wj^Kr) (15) 

1/2 

w = (?7M/7r£); n = 2-^/2 + ^1 + 4^2) (16) 
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One may evidently still use Eq. (O: 

ar = -GM/r^ - {-q^ M / Anenr) (17) 

Thus a 1/r force competes with the Newtonian one. For small M it starts to dominate at a fixed radius scale Tc, 
just as in Tohline's [H^ ^^^id Kuhn-Kruglyak's jH3| non-Newtonian gravity theories. Here Tc = 2^08 / rf . By contrast 
for M ^ Mc = ^ne/r], k « and the 1/r force scales as M^/^ and begins to dominate when the Newtonian 

acceleration drops below the fixed acceleration scale 

ao = r,V(4^Ge). (18) 

For ar ^ ao the circular velocity whose centripetal acceleration balances the 1/r force is Vc = (GaoM)^/^, precisely 
as in MOND. Thus ao here is to be identified with Milgrom's constant ao- We conclude that, with a suitable choice 
of parameters, PCG with a sextic potential recovers the main features of MOND: asymptotically flat rotation curves 
and the TFL for disk galaxies. Specifically, the choice rj = 10~^ and e = 10^'^ erg gives ao = 8.7 x lO^^cm s~^, 
Mc = 8.7 X lO^Af0 and Tc = 5.2 x 10^^ cm. Now since Tj, is larger than the Hubble scale, the Tohline-Kuhn-Kruglyak 
1/r force is comparatively unimportant. Hence for M 3> IO'^Mq we should have MOND, and for low masses almost 
Newtonian behavior. This is about right: globular star clusters at lO"* — IO^Mq show no missing mass problem. 

However, the above parameters are bad from the point of view of the solar system tests of gravity, as summarized in 
Appendix^ But the gravest problem with PCG is that it, just as AQUAL, provides insufficient light deflection [43l |. 
Here again, the conformal relation between Einstein and physical metric is to blame. TeVe S incorporates PCG's 
Lagrangian density H1Q|I in the limit of small 77 in which A becomes nondynamical. 



3. Theories with disformally related metrics 

The light deflection problem can be solved only by giving up the relation gap = e^^^gap. It was thus suggested 0| 
to replace this conformal relation by a disformal one, namely 

gc.f3 = e~^^{Agc.p + eL^Aa^,/?), (19) 

with A and B functions of the invariant g^^ ip ,1, and L a constant length unrelated, of course, to that in Eq. lIHll . 
This relation already allows to deflect light via the '4',a'4',p term in the physical metric. However, it was found j4J| 
that if one insists on causal propagation of both light and gravitational waves w.r.t. the light cones of the physical 
metric, then the sign required of B is opposite that required to enhance the light deflection coming from the metric 
gap alone. Thus the cited disformal relation between metrics, if respecting causality, will give weaker light deflection 
than would gap were it the physical metric. 

This last observation of Ref. ^3 has given rise to a folk belief that relativistic gravity theories which att emp t to 
supplant dark matter's dynamical effects necessarily reduce light deflection rather than enhancing it |34L lio. I47ll48l|. 
However, as remarked by Sanders, the mentioned problem disappears if the term ip^ai^^p is replaced by Haiip, where 
ilct is a constant 4- vector which, at least in the solar system and within galaxies, points in the time direction |5n| |. 
Specifically Sanders takes gap = e~'^^gap — 2UqU/3 sinh(2?/'). 

This "stratified" gravitation theory is reported to do well in the confrontation with the solar system tests, and 
to possess the right properties to explain the coincidence between ao of MOND and the Hubble scale But its 
vector Hq, is an a priori nondynamical element whose direction is selected in an unspecified way by the cosmological 
background. This means the theory is a preferred frame theory (although it is reported to be protected on this account 
against falsification in the solar system and other strong acceleration systems by its AQUAL behavior 0|). This is 
obviously a conceptual shortcoming which TeVeS removes, but the latter's debt to the stratified theory should be 
underlined. 



III. FUNDAMENTALS OF TeVe S 



A. Fields and actions 



TeVe S is based on three dynamical gravitational fields: an Einstein metric g^y with a well defined inverse g^'^ ^ a 

g'^^^a^p = -1, (20) 



timelike 4- vector field il^ such that 
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and a scalar field </>; there is also a nondynamical scalar field a (the acronym TeVe S recalls the theory's Tensor- 
Vector-Scalar content). The physical metric in TeVeS, just as in Sanders' stratified theory, is obtained by stretching 
the Einstein metric in the spacetime directions orthogonal to il" = g^^iip by a factor e"^'^, while shrinking it by the 
same factor in the direction parallel to iX" : 

g<.p = e-2'^(.g„0+il„%)-e2'^il„il^ (21) 

= 6^2^/3 - 2ilail/3 sinh(20) (22) 

It is easy to verify that the inverse physical metric is 

= g20^a/3 _^ 211"^'^ sinh(20) (23) 

where it" will always mean g"^iif^. 

The geometric part of the action, Sg, is formed from the Ricci tensor Rap of 5^1^ just as in GR: 

Sg - (IGttG)-! J g'^^Rapi-gf'^d^x. (24) 

Here g means the determinant of metric gap- This choice is made in order to keep TeVe S close to GR in some sense 
to be clarified below. 

In terms of two constant positive parameters, k and the action for the pair of scalar fields is taken to be of roughly 
PCG form, 

S,^-- [ [a'^h'^^c^^ac^^p + -Ge-^a^F{kGa^)] {~gf/^d^x, (25) 



2 / 2 



where /i"^ = g°''^ - il"!!'^ and F is a free dimensionless function (it is related to PCG's potential V). No overall 
coefficient is required for the kinetic term; were it included, it could be absorbed into a redefinition of a and thereby 
in k and £. Because is obviously dimensionless, the dimensions of cr^ are those of G~^. Thus /c is a dimensionless 
constant (it could be absorbed into the definition of F, but we choose to exhibit it), while £ is a constant length. 

Because no kinetic a terms appear, the "equation of motion" of a takes the form of an algebraic relation between it 
and the invariant h"'^ (f)^a<t',Pi ^^^d when this is substituted for a in Ss, the phenomenologically successful AQUAL type 
action for (f) appears. We could, of course, have written this last action directly. The present route is more suggestive 
of the possible origin of the action; for example, Ss resembles the action for a complex self-interacting scalar field 
rja ex.p{uj) / 1]) in the limit of small ?/. The term —cr^ii°'U.^(j>,a4',i3 here included in the scalar's action is new; its role is 
to eliminate superluminal propagation of the field, a recalcitrant problem in AQUAL type theories. 

The action of the vector it" is taken to have the form 

S. - J [<7"''5'^''il[a,^]%,.] ~ 2(A/i^)(5^''il^il, + 1)] {~gy/^d'x, (26) 

where antisymmetrization in a pair of indeces is indicated by surrounding them by square brackets, e.g. A^^B y^ = 
A^Bi, — AyBfj,. In Eq. H2t)|) A is a spacetime dependent Lagrange multiplier enforcing the normalization Eq. H2U|) 
(we shall calculate A later), while K is a. dimensionless constant since it" is dimensionless. Thus TeVeS has two 
dimensionless parameters, k and K, in addition to the dimensional constants G and i. The kinetic terms in Eq. H26|l are 
chosen antisymmetric not because of any desire for gauge symmetry, which is broken by the form of the physical metric 
anyway, but because this choice precludes appearance of second derivatives of IIq in the energy-momentum tensor of 
TeVe S (see next subsection) . The action is a special case of that in Jacobson and Mattingly's generalization of 
GR with a preferred frame |61j . 

In accordance with the equivalence principle, the matter action in TeVe S is obtained by transcribing the fiat 
spacetime lagrangian £{f]fj,u, 9^/", ■ • •) for fields written schematically /" as 

S„, = I C{g^,,r,r\^, ■ ■ ■){-~gY'^d^x, (27) 

where the covariant derivatives denoted by | are taken w.r.t. 5^1^. This has the effect that the spacetime delineated 
by matter dynamics has the metric g^^- The appearance of {—gY^'^ here requires us to specify its relation to (—5)^^^. 
In Appendix El we show that 

(-3)1/2 ^e-2,(_^)l/2 (28) 

By coupling to matter only through cjap, the field Uq, is totally different from the Lee- Yang 4- vector field with 
gravitation strength interaction l66|. whose existence is ruled out by the equivalence principle tests as well as by 
cosmological symmetry arguments |66l l67j . 
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B. Basic equations 

We shall obtain the basic equations by varying the total action S = Sg + Sg + Sy + Sm wth respect to the basic 
fields g°'^ , cf), a and Hq,. To this end we must be explicit about how ga/3, which enters into S',„, varies with the basic 
fields. Taking increments of Eq. H23|l we get 

g~ai3 ^ ^24,ggap _^ 2 sinh(2(/))il^(5.9^("il'') + 2 [e^'^^"'^ + 2H"'ii^ cosh(2(/.)] S(t> + 2 sinh(2(/))il("/)^(5il,, (29) 
where symmetrization in a pair of indeces is indicated by surrounding them by round brackets, e.g. A^'^B'^^ = 

1. Equations for the metric 

When varying 5* w.r.t. g"'^ we recall that 5Sg = {16iTG)^^Gaf3i—gy^^Sg°'^ {Gafj denotes the Einstein tensor of 
gap) while 

5Sm^^-fa0{-~gf'^5r'' + ■■■ (30) 

where the ellipsis denotes variations of the /" fields, and Tap stands for the physical energy-momentum tensor defined 
with the metric gap- We get 

Gap - SttG \Tap + (1 - e-4'^)il^f^(„%) + Taf^ + QaP (31) 

where 

<l>M<l>,p - ^g'^''4>,p.4',u gap - ii'^4>,f^{H(a4',p) - ^^"(f'^^dap) 
-Ge-^a^F(kGa^)gap (32) 



_ 2 
Tap = Cr 
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QaP = K (.g^''il[^.„]ii[,,0] - gap) - )^a^p (33) 

When varying g°'P in Sy we have used Eq. H20|) to drop a term proportional to gap- 

2. Scalar equation 

Variation of a in 5*5 gives the relation between a and (j)^a {F' = dF{^) / d^), 

- kGa^F - 1/2 {kGa^fF' = kfh"'^(p,a<l),p (34) 

In carrying out the variation w.r.t. <f) it must be remembered that this quantity enters in Sm exclusively through g°'^ , 
so that use must be made of Eqs. (|^ - (P?Hl : 

[<y'h''^<j,,a] .^f, = [g''^ + (1 + e-^^)il"il'^] tap (35) 

In view of Eq. H34(l this is an equation for <p only, with Tap as source. 
Suppose we define a function /x(y) by 

-f,F{f,)--fi^F'{ti)^y. (36) 

2 

SO that kGa"^ = fi{M'^h^^(j)_a(l^,f3)- We may now recast Eq. (1351) as 

[/X (kl^h^^''^,^^,.) h''U,a].j, = fcG[g"'' + (1 + e-4*)il"il'']f„^. (37) 

This equation is reminiscent of the relativistic AQUAL scalar equation [see Appendix 1X1 Eq. IjAip ]. albeit with the 
replacement g°'f^ i— > h"^ in the l.h.s. In quasistatic situations we may replace h"^ by g°'^ so that Eq. (|37|) has the 
same structure as the AQUAL equation. 
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3. Vector equation 

Variation of S w.r.t. iia and use of Eq. H29|l gives the vector equation 

Kil^"''^y^l3 + Ait" + 8nGa^il^U.i3g"'^(l)n = ^^^{1 - e-'^^)g°"'ii'^f^p (38) 

As mentioned, A here is a Lagrange muhipher. It can be solved for by contracting the previous equation with 11^. 
Substituting it back gives 

= 8TTG{l-e'^'^)[g°"'ii'^f^,p+ii"!d^ii''f^(i] (39) 

This equation has only three independent components since both sides of it are orthogonal to ila- It thus determines 
three components of il" with the fourth being determined by the normalization 1)2(1 11. Like any other partial differential 
equation, the vector equation does not by itself determine iXa uniquely. 



C. General relativity limit 

TeVeS has three parameters: k,i and K. Here we show first that in several familiar contexts the limit fc — > 0, 
i cx X cx fc of it corresponds to standard GR for any form of the function F. Many of the intermediate results 

will beuseful in Sec. Ivland I VIII We then expand on a remark by Milgrom that the GR limit actually follows under 
more general circumstances: K and £ — > oo. 

Whenever a specific matter content is needed, we shall assume the matter to be an ideal fluid. Its energy-momentum 
tensor has the familiar form 

Taf3 = pUaUp + pigap + UaUp), (40) 

where p is the proper energy density, p the pressure and Ua the 4-velocity, all three expressed in the physical metric. 
We may profitably simplify Eq. (|37|l in any case when for symmetry reasons iia is collinear with Hq,. In order that 
the velocity be normalized w.r.t. gap, we must take in that e''"da from which follows 

(jap + UaUfS = e~^'^{gap + Hq%). (41) 

Substituting this in TafS allows us to rewrite Eq. 1)37(1 as 

[fi (H2/i^>,,,0,,) = kGip + 3p) e-^*. (42) 

This form is suitable for the analysis of cosmology as well as static systems. 



1. Cosmology 

Not only important in itself, cosmology is relevant for setting boundary conditions in the study of TeVe S in the solar 
system and other localized weak gravity situations. We shall confine our remarks to Friedmann-Robertson- Walker 
(FRW) cosmologies, for which the metric can be given the form 

gapdx^dx'^ ^ -dt^ + aitfldx" + f{xf{de^ + sin^ 9 dif^)]. (43) 

Here /(x) = sinx,X,sinhx for closed, flat and open spaces, respectively. 

In applying Eq. (|37|l we shall assume that the fields 0, a and it" partake of the symmetries of the FRW spacetime. 
Thus we take these fields to depend solely on t. Also since there are no preferred spatial directions, it" must point 
in the cosmological time direction: it" = (that this is possible distinguishes il" from the Lee- Yang field which is 
ruled out in FRW cosmology Obviously this is a case where = e"^ iia', the scalar equation then takes the form 

a-^dt[a^p{-2k£'^,}>^)^] = --kG{p + 3p)e'^'l' , (44) 

2 

where an overdot signifies d/dt. The first integral is 

• . -k f* 

pi{-2kf(j>'^)(P= — G{p + 3p)e-^'l'a^dt. (45) 
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As is customary in scalar-tensor theories, we have dropped an additive integration constant; this has the effect of 
ameliorating any divergence of (/) as a ^ 0. In fact we can see that the r.h.s. of the equation behaves there as 
k{p + 3p)e^^'^t. We observe that as A: ^ with £ cx /c~'^/^, 4> will behave as k with the argument of fi staying constant. 
Thus regardless of the form of n, we have (j) ^ k. It is thus consistent to assume that </) itself is of 0{k) throughout 
cosmological history. This despite the possible divergence of 4> at the cosmological singularity, since the rate of that 
divergence is also proportional to k, as we have just seen. Recalling that kGa^ = /x, we conclude that is of 0{k'^) 
in the cosmological solutions (otherwise fi would vary with k whereas its argument stayed constant). 

Let us check whether our assumption that il" — St"' is consistent with the vector equation H38|l . The choice il" = (5t" 
makes il'"''^! = 0. For a comoving ideal fluid ii^Tap = —e^'^piia- Thus the spatial components of the vector equation 
(138(1 vanish identically, while the temporal one informs us that 

A = SvrG [a^<j}^ ~ 2p sinh(20)] . (46) 

Our previous comments make it clear that A is of 0{k). 

Turning to the gravitational equations (|31|) - (|33|l we first note that in the limit {fc ^ 0, € oc K cx fc}. Tap 

and Qa0 are both 0{k). It follows that Gap — SnGTap + 0{k). Since the difference between g^p and gap is also 
of (D{k), it is obvious that Gap = ^nGTap + 0{k) so that any cosmological model based on TeVeS differs from the 
corresponding one in GR only by terms of 0{k). In FRW cosmology TeVeS has GR as its limit when k ^ with 
£ cx fc~^/^ and K cx k. 



2. Quasistatic localized system 

We now turn to systems such as the solar system, or a neutron star, which may be thought of as quasistatic 
situations in asymptotically flat spacetime (at least up to sub-cosmological distances). We shall idealize them as 
truly static systems with time independent metrics of the form 

gap dx^'dx'^ = gtt{x^)dt^ + gij{x'') dx^'dx^ (47) 

and no energy flow. The scalar and vector equations have a variety of joint solutions. We shall single out the physical 
one by requiring the boundary condition that const, at spatial infinity, the constant being just the value of 
from the cosmological model in which our localized system is embedded. Likewise, we shall require that il" —^ dt°' so 
that the vector field matches the cosmological field at "spatial infinity" . 

We first show that il" — N^", with = the Killing vector associated with the static character of the 
spacetime, is an acceptable solution (with N = {—gap£,°'(.^)~^^^ , il" is properly normalized). Let us consider the 
expression g°'^ii^Tfj_p +ii.°'ii^iPTjp appearing in the source of the vector equation ((39|l for this choice of il". Its a = t 
component is N (T*t +iItlI*T*t^ = 0, while the a = i component is N (^g^^Tjt + il*(lI*)^Tit^ which also vanishes 

because Tjt = (no energy flow). Turn now to the l.h.s. of Eq. H39|l . Because it" has only a (time-independent) 
temporal component, il"(/)^Q. = 0, and the only nonvanishing components of ill"'''! ^^e the jt ones, and they depend 
only on the x^ . Hence ii'^'^'^^^p = so that the a ^ i components of the l.h.s. of the equation vanish. What is left 
of the a = t component is i^(il[*''^l:^ + il*iltil[*'''';/5) which vanishes by the normalization of il". Hence ii" = N 
satisfies the vector equation for any k and K. We have not succeeded in proving that this is the unique solution, but 
this seems to be a reasonable supposition. 

Now, as fc 0, the scalar equation l|57|l reduces to {iJ,h"'^(j)^a);P — 0. Multiplying this by </>(— .g)^/^, discarding 
all time derivatives, and integrating over space gives, after an integration by parts and application of the boundary 
condition at infinity, that / /i g"^(l)^a4',p{~9y^^d^^ — 0- Because for any static metric, g"^^ is positive definite and, 
when defined, /i > 0, this equation is satisfied only by </> = const, throughout. But for fc — > 0, the cosmological model 
has (/> — > 0. Hence as fc — > 0, (/> ^ in all the space. 

Returning to the full scalar equation H37|) and recalling that £ cx it is easy to see that for small but finite 

k the gradient of (j) scales as k. From the last paragraph it then follows that (f) = 0{k). These last conclusions 
are actually independent of the form of p because its argument goes to a nonzero constant in the limit fc — > 0. We 
recall [see Eq. (I34ll ] that as fc ^ 0, cr^ cx k^^. Thus the scalars' energy-momentum tensor Tap is of 0{k) (recall 
£ cx fc-3/2). From the a = t component of Eq. ^ we see that A = 0{k) + 0{K). Hence Qap = 0{k) + 0{K). 
In addition, the term in the gravitational equations (|31|l proportional to 1 — e"'^'^ is itself of 0{k)\ hence we have 
Gap — SttGTq^ + 0{k) + 0{K). Since the difference between gap and gap is of 0(0), namely 0{k), it is obvious that 
Gap — SttGTq/j + 0{k) + 0{K). Thus for quasistatic situations also, TeVeS has GR as its limit when fc ^ with 
£ cx fe~^/^ and K oc k. 
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In conclusion, the limit {fc — > 0, £ oc k~^/'^, K k} oi TeVe S is GR, both in cosmology and in quasistatic localized 
systems. 



D. Generic general relativity limit 

Milgrom (private communication) has remarked that GR actually follows from TeVe S in the more general limit 
K ^ Q and £ — > oo with k arbitrary. This is easily seen after the change of variables (/) i— > = £(/), a i-^ ct, = -y/fccr, 
whereby only gap and Ss are changed: 

~9c.i3 = e-2^*/^g„0-2il„il;3sinh(20j£) (48) 
Ss - j [ka,^h^P(j>,,a<j),j3 + -Ga,^F{Ga,-')\{-gf'-'d\, (49) 

Thus as £ ^ oo the scalar action disappears and (j)^ decouples from the theory. In addition, with if — > 0, the vector's 
action Sy disappears apart from the term with A. All this means that the r.h.s. of the Einstein equations (|31|l retains 
only the Tap and AHqH/j terms. But according to the vector equation (|38|l . from which the terms with differentiated 
4>.j, and Hq have dropped out, A ^ because (1 — e"^"^*/^) 0. Accordingly, we get the usual Einstein equations. 
Since Qafs and Qafi coincide as £ — > oo, we get exact GR. 

In this paper we shall assume that A; <C 1 and K <^ \ without restricting £. Empirical bounds on k and K are 
discussed in Sees. IIV Cl andlVl 



E. The choice of F 



Because we have no theory for the functions F{^) or ?/(/i), there is great freedom in choosing them. In this paper 
we shall adopt, as an example, the form 

plotted in Fig.^ As y ranges from to oo, ^{y) increases monotonically from to unity; for small y, /x(y) « \/yJi- 
For negative y the function /i(y) is double-valued. As y decreases from 0, one branch decreases monotonically from 
fi = 2 and tends to unity as y — > ~oo, while the second increases monotonically from ^ — 2 and diverges as y — oo. 
We adopt the second (far right) branch as the physical one. 



y 

30t 



20 - 
10 - 




-30 



FIG. 1: The function y{fj,) as relevant for quasistationary systems, < /i < 1, and for cosmology, 2 < /i < oo. 



What features of the above y{fj.) are essential for the following sections ? The denominator in Eq. (|5()() is included 
so that /i shall asymptote to unity for y ^ oo (the Newtonian limit, c.f. Sec. IIV G)l . The factor /i^ ensures that the 
MOND limit is contained in the theory (see Sec. IIV B|) . while the factor (/i — 2)^ ensures there exists a monotonically 
decreasing branch of /i(y) which covers the whole of the range y E [0, — oo) (relevant to cosmology, c.f. Sec. IVIIjl and 
only it. 
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FIG. 2: The function F{jj,) as relevant for quasistationary systems, < /i < 1, and for cosmology, 2 < /i < oo. 



Integrating Eq. H36|l with y{fi) we obtain (see Fig. 



3 M (4 + 2^-4/i2 
8 



2 ln[(l - 



(51) 



where we ignore a possible integration constant (which will, however, be useful in Sec. lVII Fl below). Obviously F < Q 
in the range /i G (0, 1) (relevant for quasistationary systems) but F > for /i > 2 (the cosmological range). Where 
negative, F contributes negative energy density in the energy momentum tensor (|32|l . Despite this there seems to be 
no collision with the requirement of positive overall energy density (see Sees. Ivl and IVII Ki . 



IV. NONRELATIVISTIC LIMIT OF TeVeS 

Sec, nil (T2I shows that in quasistatic systems TeVeS approaches GR in the limit {fc ^ 0, ^ ~ k^^^^, K ^ k}. But 
in what limit do we recover standard Newtonian gravity ? And where is MOND, which is antagonistic to Newtonian 
gravity, in all this ? This section shows that with our choice of F, both Newtonian and MOND limits emerge from 
TeVeS for small gravitational potentials, but that MOND requires in addition small gravitational fields, just as 
expected from Milgrom's original scheme. 



A. Quasistatic systems 



We are here concerned with a quasistatic, weak potential and slow motion situation, such as a galaxy or the solar 
system. As in Sec. IIII C 21 quasistatic means we can neglect time derivatives in comparison with spatial ones. Let us 
assume that the metric gaf3 is nearly flat and that |(/)| <C 1. Then linearization of Eq. H31|l in terms of the Newtonian 
potential V generated by the energy content on its r.h.s. gives gtt = — (1 + 2V^) + 0{V^). From the prescription given 
in Sec. liii(J2l il„ = -[1 + V + 0{V^)]5ta. It follows from Eq. ^ that to O(0) and 0(F), gu = -(1 + 2V + 20). 
Thus in TeVe S the total potential governing all nonrelativistic motion is $ = V" + We should remark that if 
asymptotically </> ^ 7^ 0, the gu does not there correspond to a Minkowski metric. This is remedied by rescaling 
the time (or spatial) coordinates by factors e'^" ( or e~'^'= ). With respect to the new coordinates the metric is then 
asymptotically Minkowskian. In this paper we assume throughout that |(/)c| ^ 1; Sec. IVIll shows this is consistent 
with cosmological evolution of </>. 

How is $ related to $Ar, the Newtonian gravitational potential generated by the mass density p according to 
Poisson's equation with gravitational constant G ? To relate to we first set temporal derivatives in Eq. (|42ll to 
zero which means replacing h"^(f>^a — > 9"^ 4', a'- 



kG{p + 3p)e-^'l'. 



(52) 



This equation is still exact. Next we replace g"^ — > 77"^ as well as e~^"^ — > 1. This is the nonrelativistic approximation. 
Further, to be consistent we must neglect p compared to p; keeping the former would be tantamount to accepting 
that V is not small. Thus 



Ai(H^(V0)2)v0 ^kGp. 



(53) 
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This is just the AQUAL equation Q with a suitable reinterpretation of the function ^. Now comparing Eq. H53|) with 
Poisson's equation we see that 

fc-iM|V(/)| = 0(|V$jv|) (54) 

This will be made more precise below in situations with symmetry. 

We now show that it is consistent to take V = (7$ at, with C a constant close to unity (to be determined). 
The starting point are the modified Einstein equations H31|l . With F as in H51() . F < simultaneously with 
F' < for < ^ < 1; it follows from Eq. that fi\F\ < y. Now the F term on the r.h.s. of Eq. (PU is 
—2TrG^£^^a*F{kGa'^)gai3 — ~2Trk^'^£~'^ iJ.'^F(p)gap- Similarly, since 0.i — here, the terms on the r.h.s. involving 
(p,a are of order 8TTGa'^h'^^(f>,j(j)^s gap — ^nk^'^l^'^ iiy{^)ga(3. Thus by our earlier remark the derivative terms in Tap 
dominate the F term, and by Eq. H54|l they are of order 87rfc^~^(V<I>Ar)^. But (V<i>Ar)^ is precisely the type of source 
(Newtonian gravitational energy or stress density) needed to compute the first nonlinear or 0{^n^) contributions to 
the metric. As we shall see in Sec. I VIII we need k ~ 10~^, so that if all we desire is to compute the metric to 0{^n), 
and /i is not very small, then all of Tq,^ may be neglected. 

Further, since ila = —[1 + V + 0{V^)]Sta, the il[Q,/3]^ terms in Q^p have the form (CV$Ar)^; we drop them for 
the same reason that we dropped the 0{^f^^) term in Tap- It follows that in the weak potential approximation the 
spatio-temporal and spatial-spatial components of Einstein's equations are exactly the same as in GR because the 
term proportional to 1 — e^^"^ can be dropped by virtue of the slow motion condition which suppresses the spatio- 
temporal components of Tap- The temporal-temporal component of Einstein's equations depends on A, and is thus 
another story. From Eqs. H38|l and H40|l and the observation that ii"4>,a ~ 0, 

A = Kiiaii^"'''^y,p - 167rGpsinh(2(/)). (55) 

With our Ita the first term is Kiitii^^'^/^^.p = ~KCW'^<^N+KC'^0{V<i>N'^), where by Poisson's equation V^^at = AnGp. 
Further, as we shall see in Sec. (j> is always very close to its aforementioned asymptotic value (pc (which is just (p's 
very slowly varying cosmological value). Dropping the C^OCV^n^) contribution for the same reason as above gives 

A « -8ttG[KC/2 + 2 sinh(2(/)c)]p. (56) 

Substituting this in Eq. H33|l and combining the result with the (1 — e^^'^'=) term in the Gtt equation Eq. H31|l . we see 
that (6"^"^= -I- KC/2)p replaces the source p appropriate in the weak potential approximation to GR. By linearizing 
the Gtt equation as done in GR, we conclude that 

V = {e"^''"' +KC/2)<S>N (57) 

which verifies the claim that V is proportional to <I>Ar. Indeed, since the proportionality constant here must be identical 
with C, we have C = (1 — K /2)~^e~^'^'' . Since we shall show in Sec. IVIII that it is consistent to assume |0c| ^ 1, and 
assume that if ^ 1, we shall replace C everywhere by S = 1 -I- Kl2 — 24>c. In particular 

$==S$7v + 0. (58) 

In summary, Eq. (|58l) . which is subject to corrections of 0{^n'^), quantifies the difference at the nonrelativistic 
level between TeVe S and GR, a difference which is in harmony with our conclusion in Sec. IIII G 21 We shall use 
it until we turn to post-Newtonian corrections. The condition "/i is not very small" which we imposed above to be 
able to neglect the Tap contribution to the gravitational equations is not restrictive. For the Newtonian limit we 
shall see that « 1. And when p (extreme MOND limit relevant for extragalactic phenomena), the consequent 
corrections of 0{^n^) (with large coefficient) to V are entirely ignorable because this potential is then dominated by 
4> in the expression for c.f. Eq. H59|l . 

B. The MOND limit: spherical symmetry 

First for orientation we assume a spherically symmetric situation. Then from Eq. H53|) together with Gauss' theorem 
we infer that 



In view of Eq. we have 



= (fc/47r/z)V$Ar. 

/iV<i> = V<I>Ar. 



(59) 
(60) 
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with 

/i= (S + fc/47r^)-\ (61) 

Consider the case /i < 1 for which ^i{ke^ {\V w (fc/3)i/2£| V(^| ( see Sec. IIIIE)l . Eliminating V$Ar between 
Eqs. and and defining 

we obtain a quadratic equation for /i with positive root 



fi = (fc/8^S)( -1 + ^1 + 4|V$|/ao ) (63) 

This is obviously valid only when |V$| ^ (47r/fc)^ao since otherwise /i is not small. From Eq. H61|l we now deduce 
the MOND function 



1 -l + Vl + 4|V<i>|/ao 
S l + ^l + 4|V$|/ao 

For |V$| ^ Qo (which is consistent with the above restriction since k 1) this equation gives to lowest order in K 
and (pc 

/i«|V$|/ao. (65) 

Thus if we identify our with Milgrom's constant, Eq. H6U|I with this p, coincides with the MOND formula ^ in the 
extreme low acceleration regime. Therefore, TeVe S recovers MOND's successes in regard to low surface brightness 
disk galaxies, dwarf spheroidal galaxies, and the outer regions of spiral galaxies. For all these the low acceleration 
limit of Eq. is known to summarize the phenomenology correctly. 

Now suppose |V$| varies from an order below ao up to a couple of orders above it. This respects the condition 
I V$| ^ (47r/fc)^ao. Then Eq. (|64|l shows p, to grow monotonically from about 0.1 to 0.9. Then Eq. H6UI) is essentially 
formula in the intermediate MOND regime. This regime is relevant for the disks of massive spiral galaxies well 
outside the central bulges but not quite in their outer reaches. It is known that the precise form of p makes little 
difference for the task of predicting detailed rotation curves from surface photometry. 

We see that TeVeS reproduces the MOND paradigm encapsulated in Eq. ^ for not too large values of |V<i>|/ao. 
What happens for very large | V<I>|/ao ? 



C. The Newtonian limit: spherical symmetry 

According to our choice of y(p), Eq. (|50|l . the limit /i ^ 1 corresponds to y ^ oo, that is to say |V0| oo. By 
Eqs. (|59() - (|61|l we simultaneously have |V<i>| oo and // ^ (S + k/Air)^^. Defining the Newtonian gravitational 
constant by 

Gn ^ {^ + k/47T)G, (66) 

we see from Eq. (|6U|I that V$ is obtained from V^at by just replacing G — ^ Gjv in it. In other words, in the 
nonrelativistic and arbitrarily large | V$| regime, TeVe S is equivalent to Newtonian gravity, but with a "renormalized" 
value of the gravitational constant. Now H is really a surrogate of C = (1 — K/2)^^e^'^'^''; hence for K < 2, Gn is 
positive. As mentioned, we here assume K I. 

But how close are dynamics to Newtonian for large but finite |V$|/ao ? Expanding the r.h.s. of Eq. (|50|l in the 
neighborhood of /i = 1 gives 

y = ll±+0{l-fi). (67) 

1 -/i 

We also have by Eqs. and that y = fc£^|V(/)p ~ (fc'^^^/167r^)| V^p where we have dropped corrections of 
higher order in (fc/47r). Dropping the 0(1 — /i) term in yip-) and eliminating £ in favor of ao (with S = 1) we get 

'^"^-^W ^^^^ 
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Thus to trust the approximation /i « 1 we must have |V$|/ao ^ Stt^/c ^. Using Eqs. (|68|l and H61() we obtain, 
again after dropping higher order terms in k, that 

M - 7^ U - ■ (69) 



Here the factor {G/Gm) just reflects the mentioned "renormahzation" of the gravitational constant; it is the next 
factor which interests us as a measure of departures from strict Newtonian behavior. For example, if fc = 0.03 there is 
a 5.3 X 10~^ fractional enhancement of the sun's Newtonian field at Earth's orbit where |V$| = 0.59 cm s^^. This is 
probably unobservable today. At Saturn's orbit where |V$| = 0.0065 cm s^^ the fractional correction is 4.3 x 10^^, 
corresponding to an excess acceleration 2.8 x 10^^ cm s^^ (at this point /i departs from unity by only 0.018 so that 
Eq. (|68|l is still reliable). Although this departure from Newtonian predictions seems serious, it should be remembered 
that navigational data from the Pioneer 10 and 11 spacecrafts seem to disclose a constant acceleration in excess of 
Newtonian of about 8 x 10~* cm s~^ between Uranus' orbit and the trans-Plutonian region . It is, however, unclear 
whether the correction in Eq. (jSUJ), sensitive as it is to the choice of has anything to do with the "Pioneer anomaly" . 



D. Nonspherical systems 

We now consider generically asymmetric systems. Since any system has a region where /i differs from unity and is 
variable, Eq. (|59|l is not the general solution of Eq. H53|) and must be replaced by 

V0= (fc/47r/i)(V$jv + V X /i), (70) 

where h is some regular vector field which is determined up to a gradient by the condition that the curl of the r.h.s. 
of Eq. (I7DJ) vanish. 

The freedom inherent in h allows it to be made divergenceless. Then by Gauss' theorem h must fall off faster than 
1/r^ and V x /i faster than 1/r^ at large distances. On physical grounds |V x /i| is expected to be of the same order 
as jV^Arj well inside the matter. But since the latter quantity falls off as 1/r^ well outside the matter, the curl term 
in Eq. H70|) must rapidly become negligible well outside the system. We thus expect the discussion in Sec. IIVBI to 
apply well outside any nonspherical galaxy just as it applies anywhere inside a spherical one. The interior and near 
exterior of such a galaxy, where V x is still important, must be treated by numerical methods which would be no 
different than those developed by Milgrom within the old AQUAL theory . 

Needless to say, an asymmetric system so dense that the Newtonian regime (^ approximately constant) obtains in 
its interior, e.g. an oblate globular cluster like to Centauri, can be described everywhere without an h. For in the 
interior h is not needed since even in its absence the curl of the r.h.s. of Eq. H70|l vanishes (approximately). And 
fj, begins to differ substantially from unity only well outside the system where we know from our previous argument 
that any h is becoming negligible. Hence both Newtonian and MOND regimes of the system may be described as in 
Sees. ITTBl and HVTI 

In summary, we see that the extragalactic predictions of the MOND equation are recovered from TeVe S; at the 
same time TeVe S hints at non- Newtonian behavior in the reaches of the solar system, though the effect is sensitive 
to the choice of F in the theory. 



V. THE POST-NEWTONIAN CORRECTIONS 



The upshot of the discussion at the end of Sec. lIII(T2l is that in the solar system (regarded as a static system — with 
rotation neglected — embedded in a FRW cosmological background) , Gap = 8TrGTap + 0{k)+0{K). Here we compute 
the consequent 0{k) + 0{K) corrections to the Schwarzschild metric 

g^p dx'^dx^ = _ (iz^!!!/^dt2 ^ ^ Gm/2Qf[dg' + ^^(^^2 ^ ^^^^2 ^^^2)] (^^^ 
(1 + Gm/Zg)'' 

that describes the exterior of a spherically mass m, and determine the post-Newtonian parameters of TeVe S which 
we compare with those of GR. 

Rather than just extending the Newtonian limit calculation of Sec. IIV ( ]l we start from scratch. First we write the 
spherically symmetric and static metric of the sun (inside and outside it) as 

gaf3 dx°'dx'^ = -e^dt^ + e'=[dg^ + g^idO^ + sin^ ed(p^)] (72) 
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with V = v{q) and = ^{q). Just as for metric (|71|l . outside the sun these functions should admit the expansions (ai 
and Pi are dimensionless constants) 

= l~rg/Q + a2{rg/Qf + ■■■ (73) 
= l + [3irg/Q + l32{rg/Qf + ■■■ , (74) 

where Tg is a lengthscale to be determined (see Appendix ID)) . The magnitude of the coefficient of the Vg/ g term in 
Eq. H73() has been absorbed into r^; its sign must be negative, as shown, because gravity is attractive. From the fact 
that TeVe S approaches GR for small k and K , we may infer that rg is close to 2G times the system's Newtonian 
mass. This is made precise below. 

Taking (f) — (j){g) and T^^ from Eq. H4U|I . we may write the scalar equation (|42|l as 

^-2g-(.+3.)/2[^g(.+,)/2^2^Y = fcGe-2^(p + 3p). (75) 
Here ' stands for d/dg. The first integral of Eq. (|75|l is 



Hg^ 



{p + 3p)e''^^+^'^^-^^g^dg, (76) 



where the integration constant has been chosen so that (j) is regular at g = 0. 

Supposing the matter's boundary is at g = R, we define the (positive) "scalar mass" 

m, = An r{~p + 3p)e''/2+3^/2-2^^,2^£,. (77) 

Because for a nonrelativistic fluid p ^ p, nis must be close to the Newtonian mass. In fact, as shown in Appendix IdI 
iris and an appropriately defined gravitational mass rUg differ only by a fraction of 0{Gmg j R) which amounts to 10^^ 
for the inner solar system. For g > R we may expand (f>' as 



kGr 



l + ^^+0(,-). (78) 



It is obvious from this that (j) decreases inward. Its asymptotic value, as will be explained in Sec. IVIII is positive 
and of 0{k). The decrement in (j> down to "radius" g is, according to Eq. H76|l . or its integral Eq. (|92() below, of 
0{kGms/4:iTg). In any weakly gravitating system, GrUs/g <C 1 and for strongly gravitating systems like a neutron 
star, Grus/g is still well below unity (black holes require a special discussion which we defer to another occasion). 
Thus cf) remains positive and small throughout space for all systems, and for the solar system in particular. This will 
have repercussions for the causality question examined in Sec. IVIIII 

Since we are not here interested in purely MOND corrections, we shall take p = 1 in Eq. (|78|l as well as in the 
terms in Tq^, Eq. H32|) . which explicitly involve (p derivatives. The p in the F term of Tap is not so easily disposed of 
because with our choice of and indeed with any viable one, F must be singular at /i = 1. If neglecting the F term 
in SirGTap can be justified, then using Eq. (|78|l we may compute from Eq. H32|l that for g > R 

SttGtu = 8^Gt,, = + 0{g-'). (79) 

iTT g 

Now by the approximation Hfci8|) the ratio of the F term in SnGrap to these last terms is 

8n^p'\F{pU 1287r%oV|J-(M)L . 2 
kH^GV " 3fc4|Va>A.p ^{l-^^mp)\ (80) 

which numerically does not exceed 0.04 for p > 0.99. This justifies Eq. (|79l) in any region where MOND effects 
are totally negligible. However, as pointed out in Sec. IIV CI at Saturn's orbit p already departs from unity by two 
percent. In such cases the contribution of the F term to Tq,^ must be taken into account, and its post-Newtonian 
effects compared with the MOND departure from strict Newtonian behavior calculated in Sec. IIV CI Here we shall 
only be concerned with inner solar system dynamics where p is very close to unity. Because tu is dominated by the 
derivative terms, the energy density contributed by the scalar fields is evidently positive. 
Clearly in our situation (see Sec. IIII C 2|l 

ii" = {6-^/2, 0,0,0}. (81) 
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Using this in Eqs. (|33|l and H38|) we find for g > R that 



^^!^i^±A^^!^+oi,-") (82) 



Kr 2 

With this we now turn to Einstein's equations (|31|l for all g. By virtue of 11" 's form here, the tt and gg components 
simplify to 



(t" + U'^ + 2^'/g) = SttG [(2e-4* - l)fu + r„] + Qu (85) 
<;'2 + i<;V' + + iy')/g = 8^G[f,, + r,,] + (86) 



- e 

1 „ 1 

4 



First we solve these for p > i? where Tap = 0. From Eqs. H73|l and (|74|l it follows that 

ly' = rg/g^ + {l^2a2)rgyg'' + ■■■ (87) 

^' = -Ar^/e^ + (/3i2 _ 2/32)rg 2/^3 + . . . (88) 

Substituting these together with Eqs. (|75|) . ifT^ . (f75|l and in Eqs. (|55|) - (|55|l . matching coefficients of like powers 
of 1/g, and solving the three resulting algebraic conditions gives to lowest order in k and K 

f3i = 1 (89) 

a2 = - (90) 

= - + (91) 
8 16 Snvg^ 



Using these results we show in Appendix iDlthat — 2Gmg[l+0{kGmg/ R)+0{KGmg/ R)] with uig, the gravitational 
mass, defined by Eq. ljD4|) . The relative correction here is <^ 10"^ for the inner solar system. We also remark that 
with the values l|5^ - l|nT)) the energy density contributed by Qu is positive (see Eq. (jHSl). 

For solar system tests of TeVeS we must know the physical metric 5^1,. According to Eqs. and (|HT)) . gu = 
_g20+i/^ ^ ggg/g2 = g^^/g^ svo? 9 = e~2'^+^, so we need 0. Integration of Eq. lO in hght of Eq. gives 



whereupon 



*(s) = *=-^ + 0(e-'). (92) 



The integration constant 4>c is evidently the cosmological value of 4> at the epoch in question. This value changes 
slowly over solar system timescales, so we can ignore its drift for most purposes. Thus by taking the advantage of 
the isotropic form of the metric (|72|l . and rescaling the t and g coordinates appropriately, we absorb the factors 6^'^'= 
and e~2'^'= that would otherwise appear in g^^, so that it can asymptote to Minkowskian form as expected. With this 
precaution one can calculate as if 4>c vanished. It must be stressed that this strategy works at a particular cosmological 
era. 

Accordingly 

~gu = -l + 2GArmf?-i -2/3Gw'm2f?-2 + 0(^-3) (94) 

~gg, = l + 2-iGNmg-^ + 0{g-^) (95) 

Gntu = rg/2+{kGms/AT:) (96) 

P = 1 (97) 

7 = 1 



20 



As previously, Gn is defined by Eq. (|66|) . Recalling the relations between r^, nig and m^. f Appendix ID|| . we find that 
m = mg[l + 0{kGmg/R)] + 0{KGmg/R)], i.e., in the inner solar system m and rrig differ fractionally by <C 10^^. 
Setting Vg = 2Gnig = 2Gm gives the second form of (}. Our results for (i and 7 are consistent with those obtained by 
Eiling and Jacobson [65| for the relevant case of the Jacobson-Mattingly theory. 

The (3 and 7 are the standard post-Newtonian coefhcicnts measurable by the classical tests of gravity theory ■ 
They are both unity in TeVeS, exactly as in GR (for j3 this was first noticed by Giannios). Consequently the classical 
tests (perihelion precession, light deflection and radar time delay) cannot distinguish between the two theories with 
present experimental precision. 

The /9 and 7 are not the only PPN coefficients. Future work should look at those coefficients having to do with 
preferred frame effects, as well as at the Nortvedt effect, which should not be null in TeVeS. 



VI. GRAVITATIONAL LENSING IN TeVe S 



In Sec. [3 we touched upon gravitational lensing in the Newtonian regime. Here we show that in the low acceleration 
regime, TeVe S predicts gravitational lensing of the correct magnitude to explain the observations of intergalactic 
lensing without any dark matter. First by following the essentially exact method of Ref. 0, we show this for a 
spherically symmetric structure; in nature many elliptical galaxies and galaxy clusters are well modelled as spherically 
symmetric. We then use linearized theory to give a short proof of the same result for asymmetric systems. Our 
discussion refers to lensing of both rays that pass through the system and those that skirt it, and is thus a generalization 
of the implicit result about light deflection in Sec. in more than one way. 



A. Spherically symmetric systems 



We adopt the Einstein metric H72|) : the physical metric is obtained by replacing e'^ ^^+'^4> and e*" e'^~^* in it. 
Consider a light ray which propagates in the equatorial plane of the metric (which may, of course, be chosen to suit 
any light ray). The 4- velocity i" of the ray (derivative taken with respect to some suitable parameter) must satisfy 



(99) 

From the metric's stationarity follows the conservation law e'^+^'^i — E where E is a, constant characteristic of the 
ray. From spherical symmetry it follows that e''~'^'^g^(p = L where L is another constant property of the ray. Let us 
write g = {dg/d(p)ip. Now eliminating t and ip from Eq. I|99|l in favor of E and L, and dividing by E^ yields 



,-1.-20 



+ {b/gye-^+'n9-'idg/dipy + l]=0, 



(100) 



where b = L/E. By going to infinity where the metric factors approach unity one sees that b is just the ray's impact 
parameter with respect to the matter distribution's center at p = 0. This last equation has the quadrature 



if : 



(I)' 



1 



-1/2 dg, 
Q ' 



(101) 



Were the physical metric exactly flat, this relation would describe a line with (p varying from to tt as p decreased 
from infinity to its value gtum at the turning point, and then returned to infinity. Hence the deflection of the ray due 
to gravity is 



- 1 



^/^dg 



Q 



(102) 



This last integral is difficult. So let us take advantage of the weakness of extragalactic fields which allow us to 
assume that ^ and (j) are all small compared to unity. Then the above result is closely approximated by 



da 



(1 + 5 



1/2 dg, 

Q 



(103) 



The rewriting in terms of an a derivative allows us to Taylor expand the radical in the small quantity (, ~ v — Akp 
without incurring a divergence of the integral at its lower limit. The zeroth order of the expansion yields a well known 
integral which cancels the tt. Thus, to first order in small quantities 



2d_ 

b da 



A(f>)gdg 



(104) 
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At this point it pays to integrate by parts: 
2 d r 



bda 



Q^OO 



lim ^){g^ ~ ab^f'^ ~ ~ v' - 40')(e' - ah^f'^dg (105) 



a = l 



Since v, ? and 4> all decrease asymptotically as , the integrated term, being a independent, contributes nothing. 
Carrying out the a derivative, and introducing the usual Cartesian x coordinate along the initial ray hy x = ±{0^ — 
6^)^/2, we have 

A,^V°° ^'-^'^'^'dx. (106) 

A factor 1/2 appears because we have included the integral in Eq. (|105|) twice, once with g decreasing to, and once 
with g increasing from h. The integral is now performed over an infinite straight line following the original ray. 

The difference between GR with dark matter and TeVe S in this respect is that with dark matter one would have 
= and would compute u and ^ from Einstein's equations including dark matter as source, whereas in TeVe S one 
has a nontrivial 0, and computes u and on the basis of the visible matter alone. 

We may simplify the above result by means of Einstein's equation l|86|l . We shall neglect the and q'v' terms 
because they are of second order, and thus smaller than v' / g by a factor G ■ mass/g which amounts to , with v the 
typical orbital velocity in the system. Using the residual terms we eliminate <j' from Eq. (|106|l : 

v' + 26' 1"°° ~ 

Aip = b ——I- dx - AnGb / (Tgg + + eee/SvrG) dx. (107) 

J —OQ Q J —OC 

Now by Sec. IIV Al ly ^ 2V + 0(V'^) and ^ = V + <j). Hence with fractional corrections of 0{V'^), 

/OQ P '^^^ 

— dX-4TTGb / {fgg + Tgg + egg/8TTG)dX. (l08) 

-00 Q J —oo 

The first integral here depends exclusively on the potential $ which determines nonrelativistic motion. That is, the 
observed stellar or galactic dynamics will uniquely fix this part of Atp. For this reason the first term makes the same 
predictions for lensing by nonrelativistic systems in TeVeS as in GR (where ^ — the last calculated assuming 
dark matter). We next show that for nonrelativistic systems the second integral is negligible. 

In astrophysical matter the radial pressure Tgg is of order p times the local squared random velocity of the matter 
particles (stars, gas clouds, galaxies). Thus / Tggdx = {v"^) J pdx with (w^) a suitably averaged . But by Poisson's 
equation AnGp — V • V<i>Ar ^ ^n' / g — P'^' / g where we have also used Eq. jnOJ. Thus the term with the integral 
over Tgg is smaller than the first term in Eq. p()8|l by a factor of 0{fl{v'^)). In GR (for which effectively p. = I) this 
factor is no larger than 10^^ for all extragalactic systems which have a missing mass problem; in TeVe S it is even 
smaller because typically /i <C 1 for such systems. 

Turning now to Tgg we recall from Sec. IIV .^ that in the quasistatic situation in question, the F part is dominated 
by the term quadratic in derivatives. Using Eqs. (|59|l - Ht)0|l we work out that AirGrgg « {kjl/STrp)^'^^' ■ Evidently 
(j)' ^ ^/g, and since <I> = 0{v'^) and {kp/Sirp,) < i, the contribution of Tgg to the second term of Eq. H108|) is no 
larger than that coming from Tgg. 

Finally we note that the A term in Qgg vanishes in a quasistatic situation because then Hq ~ —(1 + ^N)^at- And 
from this last formula we estimate \Qgg\ ~ ^K{^]y')^ Kjl^l^^'l/ g. Since p < 1 and by Sec. El we must take 
K < 10~2, it is clear that the contribution of Qgg is much smaller than that coming from Tgg. From all the above the 
light ray deflection in TeVe S is 



A(^ = 26[l + 0(/V)] / —dx. 

J —OQ Q 



(109) 



In GR with dark matter the same formula is valid with 0{pv'^) replaced by 0{v'^). Since these corrections are beyond 
foreseeable accuracy of extragalactic astronomy, it is clear that for given dynamics (given <&), both theories predict 
identical lensing. We shall elaborate on this statement shortly. 



B. Asymmetric systems 



We now turn to systems with no particular spatial symmetry. The weakness of the gravitational potentials typical 
of nonrelativistic systems entitles us to use linearized theory |6J| in which the metric is viewed as a perturbed Lorentz 
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metric: 

gaf3 = VaP + haf3 - - rjaP V'^^hjS (110) 

2 

with \hap\ ^ 1. By small coordinate transformations one enforces the gauge conditions rj^^h^s^p = 0; as a consequence 
to first order in the h fields 

Ga0 = --Ti^^d^d5Kp. (Ill) 

2 

SO that Einstein's equations take the form of wave equations in flat spacetime with the r.h.s. of Eq. (|31|l as sources. 
Of course there are motions and changes in galaxies and clusters of galaxies, but the associated changes in the metric 
are mostly very slow. Thus we confine ourselves to quasistationary situations where we can drop time derivatives (but 
not yet the gu components since galaxies do rotate). This tells us that 



Qtt- (112) 



Gu = — htt = SttG ftt + 2(1 - e-4'^)il^f^tit( + m 

2 L J 

The various parts of the source here were explored in Sec. IIV Al from that discussion it follows that 

htt = -4V = -4S$w. (113) 

In regard to the spatio-temporal source components of Eq. H31|l . we observe that the Ta is an 0{v) below Tu 
(momentum density is velocity times mass density). Further, the dominant contributions to th are hu multiplied by 
(j'^r]^^(j) j(j) f. and by {GjP'^u^F. Of these the first dominates (see Sec. IIV A^ . and it is small on the scale of p both 
because it is of second order (c.f. Sec. 0), and because \^ti\ ^ 1- We can guess that il.; is at most of order hu (it 
would vanish in a truly static situation), and since by Eq. I|56|) A is below SttG/S by factors of 0{K) and ©(^c), we 
see that the Ailjili term contribution to is small compared to SirGp. Similarly, the Kg'^'^U[fj_ t]^[^ i-^ contribution 
to Qti, being of second order in Vi and first order in hu, or first order in Vi and first order in hti,j (aside of carrying 
the small coefficient K), must be very small. We conclude that the source of the spatio-temporal Einstein equation 
can be neglected, so that to the accuracy of Eq. hu ^ 0. 

Things are similar for the spatial-spatial components. We have already remarked that Tij is an 0{v'^) below Tu- 
The Tij consists of a term quadratic in (jj^i and one with a F factor which has been argued to be smaller. Hence 
is small. Again the iiT^^'^ilj^, contributions to 0^^ are quadratic in Vj and suppressed by the K coefficient, so 
they are also small. And the A, which we remarked above to be small, is multiplied by two factors ha, and so is also 
small. So by the same logic as above we neglect the sources of the spatial-spatial components hij and conclude that 
hij ~ 0. 

Substituting all these results in Eq. (|110|l we obtain 

gap = (1 - 2V)T]ap - iVSatSpt. (114) 

The absence of gu in this approximation makes the situation truly static (rather than just stationary); hence il" — (5f . 
Calculating the physical metric from Eq. (|22|) with e*^"^ w 1 ± 20 we have 

gap = (1 - 2F - 2(j))r,ap - 4(y + ^)SatSpt (115) 

which is equivalent to 

gapdx"dx>^ = -(1 + 2$)di2 + (i _ 2^)S^jdx^dx^ (116) 

with — V + d> StS in 

Metric (|116|) has the same form as the GR metric for weak gravity |6^. Thus in TeVeS just as in GR the same 
potential governs dynamics and gravitational lensing. This accords with the conclusion of Sec. I VI Al for the spherically 
symmetry case. What does this mean in practice ? In GR <I>'s role is played by the Newtonian potential due to the 
visible matter together with the putative dark matter; in TeVe S is the sum of the scalar field and the renormalized 
Newtonian potential generated by the visible matter alone. These two prescriptions for $ need not agree a priori, 
but as we argued in Sec. lIVBl nonrelativistic dynamics in TeVeS are approximately of MOND form, and MOND's 
predictions have been found to agree with much of galaxy dynamics phenomenology. We thus expect TeVeS's 
predictions for gravitational lensing by galaxies and some clusters of galaxies to be as good as those of dark halo 
models within GR. But, of course, the early MOND formula (^3), and TeVeS with our choice l|51|) for F(/i) both claim 
that asymptotically the potential $ of an isolated galaxy grows logarithmically with distance indefinitely. Dark halo 
models do not. So TeVe S for a specific choice of F is in principle falsifiable. Dark matter is less falsifiable because 
of the essentially unlimited choice of halo models and choices of their free parameters. One should also remember 
that gravitational lensing affords the opportunity to map the $ to greater distances than can dynamics; for unlike 
the latter, lensing can be measured outside the gas or galaxy distribution. Using this $ both GR and TeVe S would 
predict the same dynamics for stars or galaxies, while disagreeing on the implied distribution of mass. 
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VII. COSMOLOGICAL EVOLUTION OF </> 
A. Persistence of cosmological expansion 

This section (where we write (p rather than 0c) shows that for a range of initial conditions, FRW cosmological 
models with flat spaces in TeVeS expand forever, have < </> <C 1 throughout, and their law of expansion is very 
similar to that in GR. The second point is crucial for our discussion of causality in Sec. IVIIII 

First using Eq. (|22|l we transform metric H43|l to the physical metric 

go,pdx°'dx'^ = -dp +h{if[dx^ + f{xf{de^ +sm^edif^)\, (117) 
di = e'^dt; a = e~'^a. (118) 

In what follows we take the initial moment, conventionally written as i = 0, at the end of the quantum era with a(0) 
a very small scale; furthermore we take the zero of t to coincide with t = Q. For illustration we assume the initial 
conditions ^(0) = (an overdot always denotes d/dt) and < <j)Q = (/)(0) <C 1. Hence a also starts off from a very 
small scale, ap, and can only increase initially. 

We now show that the spatially flat (/(x) = x) FRW models in TeVeS persist and cannot recoUapse, i.e. d has no 
finite maximum. As in Scc. lIIICll we have il" — St" which causes H^"'^^ to vanish. As a consequence Qap — ~A(5^(5^ 

with A given by Eq. Since (j> = (j>{t), Eq. ^ gives m = 2a'^^'^ + G{4:e'^)-^a'^F(p). As mentioned in Sec. lIIICll 

iiP fafi = -pe^'f'Ha. Using g'^'^iiciip = -1 gives us fu + (1 - e-'''^)il"f'„((ilt) = {2e-'^'t' ~ l)pe^'t'. Substituting all the 
above in the tt component of Eq. H31I) . we get the following analog of Friedmann's equation: 

With the choice H5()|l for y(/i) we have p > 0, y{lJ.) < and F > in the cosmological domain. Thus the scalar fields 
contribute positive energy density and the r.h.s. of Eq. (|119|l is positive definite (/5 < is physically unacceptable). 
It follows that d cannot vanish for any t, so that by our earlier remark it must always be positive. Now the relations 
(|118|l imply that 

dd/dt^ e-^'''{d~a(j)). (120) 

We shall show in the sequel that although </> can be positive, it is always the case that |0| <C d/a. As a consequence 
dd/di is always strictly positive: in TeVeS a FRW model with flat spaces cannot recoUapse. 

The fact that is given by an integral over time [see Eq. H45|l ] means that in a cosmological phase transition, where 
p may change suddenly, (p (and of course (p) will nevertheless evolve continuously in time. It follows that F will also 
evolve continuously in time [see Eq. (|36|) ]. A consequence of Eq. (|119|l is that any jump in p will be reflected in a 
similar jump in (a/a)^ or in the square of the Hubble function H = dd/dt. 

B. The proto-radiation era 

Contemporary cosmology regards the inflationary era as preceded by a brief radiation dominated era, the proto- 
radiation era, in which the physical scale factor d expands by just a few orders following the quantum gravity regime. 
As in any radiation dominated regime, here the equation of state is j5 — 3p with both p and p varying as d~^. It 
follows from Eq. (|45|l that throughout the era 



k 



= I G~pe~^'^a^dt, (121) 







Because in the cosmological regime /i > 2, we have </> < throughout this era. Thus as promised dd/dt in Eq. H12()(l 
is positive. Using the constancy of (Gp)^/^a^e~^'^ we can now write 

,0^-MM^r(Gp)V^a.. (122) 
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Tentatively assuming that jt^l <C 1 throughout the era we may, according to Eq. H119() . bound both instances of 
(Gp)^/^ from above by (3/87r)^/^a/a. The consequent integral is then trivial, and since aq is essentially zero we get 

< (3fc/87r) (a/a). (123) 

Thus \4>\ < (3fc/167r)(d/a); since fc < 1, we have by Eq. I|12U|1 that da/di « d. 

We can now show that the cosmological evolution during the proto-radiation era is very similar to that within GR. 
For the choice H51|l both F and F' are positive in the cosmological domain (see Fig. ^ . It follows from Eq. H36|l that 
fj.'^F < —fiy (recall that y < 0), so the last term on the r.h.s. of the Friedmann equation is less than half the second. 
Next we use y = — 2/c£^(^^ to infer from Eq. (|123() that 



But this means that the scalar field contributions to the Friedmann equation are small compared to its l.h.s. Specifi- 
cally, to within a fractional correction of 0{k/16) (actually smaller than this because /i will turn out to be large), the 
relation between H and p is the same as in GR. 

The fact that the scalar field contributions to the Friedmann equation are small compared to its l.h.s. also means 
that inequality H123|) is nearly saturated, as must be its kin (|124|) . Then 

p^\yip)\^-{3k/A7r)\d/afao-'. (125) 

6 

But a/d is a very short scale (in standard cosmological models 10~^^ s in the proto-radiation era) while 

ao~^ 3 X 10^^ s. Thus ii^y{p) ^ 1. Since by Eq. H5U|) this is possible only for /i ^ 1, we can sharpen our earlier 
conclusion from Eq. H123|) : \4>\ ^ (3fc/87r)d/a. Now it is even clearer that a and d (as well as t and i) are essentially 
equal, so that the expansion in this era proceeds just as in GR. Further, integrating this last inequality gives 

\(i3pr - 0o| < (3fc/87r) ln(ap,/ao), (126) 

where the subscript "pr" stands for the end of the proto-radiation era. Since this era spans just a few e-foldings of 
the scale a, the logarithm here is of order unity. Hence </> is almost frozen at its initial value (/)o, provided this last is 
not extremely small. By choosing as initial condition Q < (j)^ 1, as we proposed, but avoiding extremely small ^o, 
we get < ^ 1 throughout the proto-radiation era, as assumed earlier. Thus our assumption was consistent. 

C. The inflationary era 

The equation of state during inflation is p = —p = const. Then (|45|) tells us that 

^0=A^ Gpe-^^a'^dt + pj,r^pr[^)^ ■ (127) 

The integration constant prefacing the last term is fixed by the condition that p and </> be continuous through the 
proto-radiation inflation divide. It is clear that after rapid expansion has suppressed the last (negative) term here, 
(j) becomes positive. Because p is constant, we may pull a factor {GpY^^ out of the integral. Then by Eq. H119|l and 
assuming everywhere that e~"^ « 1 (which we verify below), we have (G/5)^/^e~^''^ < {S/SirY^^d/a both in and outside 
the integral. Thus 

3/cd /■* 2 • /'0-pr\^ 

p4> < - — J / a ddt + ppr(ppr[ (128) 



where we have used Eq. (|123|l as an equality as the end of the proto-radiation era. Thus during inflation 

- {3k/8TT){d/a)pr < p<i)< {k/8n){d/a). (130) 

The l.h.s. here comes from the last term in Eq. (|127|l in light of inequality (|123|l . In the passage from the proto- 
radiation era, which involves a phase transition, p can change by a factor of order unity, but then settles down to a 
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constant. Thus by Eq. H119(l aja remains at least of the same order of magnitude as {a/a)pr. Hence inequality H130|) 
translates into one of the same form as 1)1231) but valid during inflation. As in Sec. lVllEl this tells us that da/dtK, a 
also during inflation. And the argument following inequality H123|l can now be repeated to show that the — /iy and 
fi^F terms in Friedmann's equation amount to relative corrections of 0{k/16) (actually smaller), so that inflation in 
TeVe S proceeds very much like in GR. 

Repeating the argument leading to Eq. (|129|l in light of this last conclusion and the added realization that the 
terms disappear very rapidly, we conclude that during the > part of inflation, inequality (|123|l is very nearly 
saturated. One can then rederive Eq. 1)1251) as in Sec. IVllljl Because the inflation timescale is again very short 
compared to ao^^, the argument yielding Eq. (|126|l can be repeated with slight modifications to show that during 
inflation /i ^ 1, and consequently 

4>i - (j}pr <^ i3k/87:)ln{ai/apr), (131) 

where a subscript "i" stands for the end of inflation. Thus, although in standard models inflation can span up to 
70 e-foldings of a, the r.h.s. of this inequality is very small compared to unity. We conclude that inflation manages 
to raise (j) above its value at the end of the proto-radiation era by a very small fraction of unity. This justifies our 
replacement of e"*^ by unity in deriving Eq. H129|) . 

In what follows we shall denote by Hi, fii, 4>i and the values of the Hubble parameter, /i(— 2/c£^(/)^), (j) and (j>, 
respectively, at the end of inflation, t = ti, where a — ai. 

D. The radiation era 

In the ensuing radiation era the equation of state switches back to "ip — p with both p and p varying as a~^. Thus 
the integral in Eq. H45|l is 

Gpe-2Vd^ + ^,0,(^)^ (132) 

with the integration constant Pi4>i set so p4) at the radiation's era outset equals that at inflation's end. Although 
initially (/) > 0, clearly the integral will eventually dominate the last term making cf) negative thereafter. 

Now according to Eq. (|129() . < {k/87T){a/a)i. Due to the approximate continuity of (a/a) across the inflation- 
radiation eras divide [which itself follows from the approximate continuity of p and Eq. H119() ]. and from the fact that 
(d/a) falls off no faster than (ai/a)^ in the radiation era, Eq. 1)132)) gives 

< {k/8TT){d/a)i {ai/af < {k/8TT){a/ a). (133) 

On the other hand, from pa^ — const, we can move a factor (G/3)^/^a^e~^'^ out of the integral in Eq. 1)132)) . Using 
again (GpY^^ < (3/87r)-'^/^d/a from Eq. ()119)) (if we assume provisionally that e^"^ « 1) both in and outside the 
integral, we have 

Ai</'>-(|^)^ (<i/a)adi + /i,0,(^)'. (134) 

The integral is a{t) — a^. Hence 

p,:j) > (-3fc/87r)(l - a,/a)(d/a) + ^i(/),(a,/a)^ > -(3A:/87r)(d/a) (135) 

In view of Eqs. 1)133)) and 1)135)) . inequality ()123)l is again valid here. Because p > 2 we get again from Eq. ()12())l that 
da/ da « d/a. We may now reproduce inequality ()124)) and show as in Sec. lVIIBl that to within a fractional correction 
of 0(^/16), the relation between H and p is the same as in GR. 

Because of this last result, Eq. 1)133)) and the rapid decay of a^/a in Eq. I)135)) . we may conclude that when < 0, 
inequality 1)123)) is nearly saturated. We may then rederive Eq. 1)125)) as before. Now in conventional cosmology at 
redshift z during the radiation era H 3x 10~^''(1 + z)^ s~^, which by previous inference closely approximates d/a in 
our model. We thus obtain /.t^|y(/i)| « 5 x 10^^A;^(1 + z)"*. Takiner k ~ 0.03 on the basis of Sec. lIV Cl we see that at the 
end of the radiation era {z « 10^), p'^\y{p)\ « 4 x 10^ which corresponds to /Lt « 10. For earlier times /i cx (1 + z)^^^ 
so that it rises to 10^^ at the beginning of the era at z « 10^^. Going back to inequality ()123)) we see that in the 
last three e-foldings of the era (j){t) — > —8 x 10"'' with the previous 50 e-foldings contributing an even smaller 
decrease. Our assumption that e~'^ « 1 was evidently justified if (j>o is taken small compared to unity, yet sufficiently 
positive to keep (j){t) positive throughout the era. 

We shall denote by /i^, 4>r and 0^ the values of p{~2ki'^(j)^), (p and 0, respectively, at the end of the radiation era, 
t = tr where a = a^. 
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E. The matter era 

In the matter era p « and p varies as 5,^'^. Integrating Eq. 145(1 gives, c.f. Eq. H132|l 

p^ = ^ f Gpe^^'l'a^dt + Hr^r(—Y ■ (136) 

It is clear that (j) continues to be negative throughout the matter era. Using pa^ e~^'^ = const, and setting henceforth 
— 1 (whose consistency will be checked below), we explicitly perform the integral in Eq. H136|) from U to t: 

Il4' — kGp{t — tr) + Pr4>r{lr/o-)^ ■ (137) 

Integrating the inequality (Gpa^)^/'^ < (3/877)1/20^^0 coming from Eq. l(TT^ we get {Gp)'^/'^{t-tr) < (2/3)(3/87r)i/2. 
Both together give Gp{t — tr) < (d/47ra), which when substituted in Eq. ((137(1 finally gives 

p(f) > (— fc/87r)(d/a) + Pr'Pr {ar/af. (138) 

Now according to Eq. ((135(1 /ir^r > (— 3fc/87r)(d/a)j.. Thus at the beginning of the matter era, where a — ar, the 
lower bound on the second term on the r.h.s. of inequality 1(138(1 maybe as much as three times larger in magnitude 
than the first term, yet it decays as a~'^ while the first term cannot do so faster than a~^^^ [see Friedmann's equation 
((119(l ]. Hence within about one e-folding of a, the first term comes to dominate the r.h.s., and over most of the matter 
era 

< (fc/87r)(d/a). (139) 

From this follows a tighter version of bound p24l) which again demonstrates that the scalar field terms in Einstein's 
equations are rather negligible. The fact that 11139(1 may be exceeded by a factor of a few early in the matter era is 
no reason to exclude that epoch from the just mentioned conclusion: the rather large p at the end of the radiation 
era {p ~ 10 ) — and a bit beyond — acts to suppress that factor. Using by now well worn logic we conclude that in the 
matter era as well, the relation between H and p is almost the same as in GR. 

Integrating inequality 1(139(1 with the use of /i > 2 (the first e-folding's relatively larger contribution is suppressed 
by the larger p which holds sway then), we get 

(j){t) -(t)r> -(k/ien) ln(a/a^). (140) 

Because the matter era thus far has spanned nine e-foldings, has decreased by less than 0.0054 during this era. 

Note that we have not addressed the cosmological matter problem. In TeVe S the expansion is driven by just p, 
the visible matter's density, whereas the observations require that the source of Friedmann's equation which falls 
off like should be larger by a factor of perhaps 6. There are at least two possible avenues for dealing with 
this embarrassment. First, we have stuck to a particular F(p); possibly a more realistic F(p) would change late 
cosmological evolution enough to resolve the problem. Second, we have insisted on being small. This is a consistent 
solution as we have shown, but it is perhaps not the unique solution. Plainly nonegligible values of cj) can affect the 
Friedmann equation significantly. 

F. The accelerating expansion 

Lately data from distant supernovae indicate that in recent times {z < 0.5) the cosmological expansion has began 
to accelerate, namely, that dH/dt > —H^. The data are best interpreted in GR by accepting the existence a positive 
cosmological constant A « '^Hf^^^y [g^ . One can incorporate such accelerating epoch in the TeVe S Einstein equations 
(PT|l by adding to p'^F{p) — purely phenomenologically — a constant (/^-independent) term of magnitude « Ki^f^jl-n. 
Such constant part, which corresponds to the integration constant involved in solving Eq. (|S^ . leaves y{p) unchanged, 
merely shifting the curve for F(^p) in Fig. [3 up. Furthermore, according to Eq. I(62() and the empirical connection 
ao ^ Hq 0, the added constant is ~ Sk^/ldn^, that is very small. It cannot thus affect the discussion in earlier 
sections, and in particular F continues to make a positive contribution to the energy both in static systems, and in 
cosmology. 

The appearance of the cosmological constant in F has almost no effect on the value of cj). To see why note that 
A does not directly affect the scalar equation ((42(1 . but only the Friedmann equation ((119(1 . Hence Eq. ((137(1 is still 
valid. As the expansion accelerates, a begins to grow exponentially with t. Both terms on the r.h.s. of Eq. 1(137(1 thus 
fall off drastically, and (j) becomes "stuck" at the value it had soon after the onset of acceleration. Consolidating the 
results of Sees. IVII BIVII El with our conclusion we see that the range of initial conditions 0.007 < 0o ^ 1 insures 
that (j) > and e"^ w 1 throughout cosmological evolution. 
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VIII. CAUSALITY IN TeVe S 

TeVeS's predecessors, AQUAL and PCG, permitted superluminal propagation of scalar waves on a static back- 
ground. In the case of PCG with a convex potential this occurs hand in hand with an instability of the background, 
so it is unclear if true causality violation occurs. How does TeVe S handle causality issues ? 

The question is complicated here by of the existence of two metrics, gap and gap^ whose null cones do not coincide 
(except where (/) = 0). Which of the two cones is the relevant one for causal considerations ? We shall take the 
view that since common rods and clocks are material systems with negligible self-gravity, the coordinates to which 
the Lorentz transformations of special relativity refer are those of local orthonormal frames of the physical metric 
(jafj and not of gap- It is by ascertaining that in no such physical Lorentz frame can physical signals travel back 
in time that we shall certify the causal behavior of TeVeS. Now Lorentz transformations involve a parameter, the 
critical speed "c". We shall identify this with the speed of electromagnetic disturbances so that, as customary, the 
speed of light is the same in all Lorentz frames. Since we have built special relativity into TeVe S by insisting that all 
nongravitational physics equations (including Maxwell's equation) take their standard form when written with gap^ 
this procedure is consistent. In fact, all signals associated with particles of all sorts are subluminal or travel at light's 
speed with respect to cjap- 

There remains the question of whether gravitational perturbations (tensor, vector or scalar) can ever exit ga/3's 
null cone. The analysis given below is quite different for tensor and vector perturbations on the one hand, and scalar 
perturbations on the other. One point in common, however, is that causality is guaranteed only in spacetime regions 
for which > 0. As shown in Scc. lVIII there is gamut of reasonable cosmological models for which <p is indeed positive 
throughout the expansion. 

A. Propagation of tensor and vector disturbances is causal 

The characteristics of both Einstein's equations (|31|1 and the vector equation (|38|l lie on the null cone of gap because 
all terms in them with two derivatives are the usual ones in Einstein's and gauge field's equations. Accordingly, we 
do not expect metric and vector perturbations to travel outside the null cone of the Einstein metric gap- However, 
the interesting question is rather what is the speed of a wave of this class in terms of the physical metric gap ? 

In the eikonal approximation the wavevector of metric perturbations, that is the 4-gradient of the characteristic 
function, will satisfy g^^Hai^p — 0. Hence Eq. gives 

g'^^Kai^p - 2(il"K„)^ smh{2(j)) = 0. (141) 

We consider a generic situation where il" may have both temporal and spatial components. The normalization l|2()(l 
implies by Eq. H22I) that cjap'd"'^^ — — e^"^. Thus in an appropriately oriented local Lorentz frame, £, of the metric 
cjap we may parametrize il" by 

il" = e'^(l - F2)-V2{i^_V",0,0} (142) 

with —1 < y < 1. This V is actually the ordinary velocity (measured by the physical metric) of £ w.r.t. the 
privileged frame in which the matter is at rest (whether in cosmology or in a local static configuration) , namely that 
in which il" = {e*^, 0, 0, 0}. This is evident by considering a Lorentz transformation from the matter rest frame to the 
coordinates appropriate to frame £. 

In view of the above, Eq. H141() reduces to 



= Au'^ + 2Bn\\uj + Dk\\'^ - {1-V'^)ka_'^ (143) 

A = eZ-^-y^ (144) 

B = V{e^'l'-l) (145) 

D = -l + T/2g40 (--^^g) 



with to = —Kt and /ty and kj_ the spatial components of Ka collinear and normal to iii (the space part of Hq), 
respectively. For arbitrary V (|143|l is an anisotropic inhomogeneous dispersion relation (uj depends on position 
through (f) as well as on direction of the wavevector). However, in the rest frame of the matter {V — 0), it is isotropic 
(though still position dependent through 0) with group (or phase) speed equal to 



^^0 = e-^t 



(147) 
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The condition for tensor and vector perturbations not to propagate superluminally {vq < 1 as judged in the physical 
metric) is thus that </> > 0, which as we saw, is satisfied in a wide range of cosmological models (see Sec. lVll() as well as 
quasistatic systems embedded in them (Secl^. Normally this conclusion could be carried over to all Lorentz frames 
without further calculations. But because TeVe S admits a locally privileged frame, that in which 11" = e'f'{l, 0, 0, 0}, 
we investigate this conclusion in more detail for any V'^ < 1. 

Solving Eq. Ij 143(1 for lu gives 

uj ^ {-Bk\\ ± S)A-^, (148) 
S = {Ck\+A{1-V^)k\Y''', (149) 
C = - AD = [l-V'^fe'^'^. (150) 

The condition > makes A here strictly positive. It is possible for the above expression for lo to change sign, so for 
given K we must agree to always choose the branch of the square root that makes u) positive (negative to with opposite 
sign K is the same mode, of course). In what follows we call the modes with upper (lower) signs of the square root 
+ (— ) modes. For the components of group velocity collinear and orthogonal to ili, respectively, we derive 

v|| = dijj/dK\\ = {-B±CS-^K\\)A-^, (151) 
v_L = dLo/dK,^=±{l~V'^)S-^n^. (152) 

Since these expressions are homogeneous of degree zero in /«, there is no dispersion, but for V ^ Q the propagation is 
anisotropic. For small </) one has analytically 

u = 1 - 2(1 ± l/cos7?)2 (1 - V^)-^(f) + 0{(t)^) (153) 

where v = (v| + |v_l p)!/^ and d is the angle between k and it;. Thus for moderate V the group speed v is subluminal, 
but obviously formula l|153|l becomes unreliable for V close to unity. 

For arbitrary V it is profitable, as remarked by Milgrom, to write v in terms of lo. In fact a straightforward 
calculation gives 

l-v^ = S-^C{kI+ Ki-uj^), (154) 

from which it is clear that v can become supcrluminal only if the (isotropic) phase speed ti-'(K| + m"]^)^^/^ docs the 
same simultaneously. Since the latter was found subluminal at y = 0, we have only to ask if there is some V < \ for 
which w = (k| + k\Y/'^] we might then suspect there is supcrluminal propagation for larger V . Suppose we substitute 
this last value of lo in Eq. (|143|l together with those of A, B and D. Collecting terms one can put the condition for 
the transition to superluminality in the form 

(e^^-l)(l/^ll + yi^fT^)' = 0. (155) 

As we saw in Sec. lVIII for a broad class of cosmological models > throughout the expansion, and as Sec.lVltestifies. 
variation of in the vicinity of localized masses embedded in such a cosmology is far short of what is required to 
turn the sign of (j>. It is thus clear that even in the case < 0, condition (|155|l cannot be satisfied for V < \. Hence 
supcrluminal propagation of vector and tensor perturbations is forbidden. 

How does V vary with V 1 When k±_ ^ 0, we find numerically the following behavior. For the + mode with < 0, 
vy < for all V ^ and after experiencing a shallow maximum at modest V , v reaches a minimum at V very near unity, 
which is the deeper and farther from V = I the larger |/«_L//«y|. As V grows further, v rises and approaches unity for 

V ^ 1. If ^ii > 0, V|| starts positive for small V but eventually turns negative at a rather large V which grows with 
|Kj_/Ky|. As V grows further, v reaches a minimum, which gets shallower with growing |Kj^//«y|, and then begins to 
rise. At a critical V the positive Ky + mode terminates. However, the — mode with negative k\\ takes over onward 
from the critical V; it features vy < 0, and for it v rises with V and approaches unity as — > 1. The — mode with 
Ky > is always unphysical. 

For n± = and Ky < the + mode has vy < throughout, and v rises monotonically with V approaching unity as 

V 1. For /ty > that same mode has vy > and v decreasing with increasing y up to a = Vc ~ e^^"^ at which 
point both vy and v vanish. The terminated sequence is continued by the — mode with Ky < for which vy < and 

V rises monotonically with V from zero at V = Vc and approaches unity as V 1. 
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B. Propagation of scalar perturbations is also causal 

The terms with second derivatives in the scalar equation H37() have a nonstandard form reminiscent of those in 
relativistic AQUAL (see Appendix^. Do scalar perturbations propagate across ga/j's null cone, that is do they 
travel faster than electromagnetic waves ? We now show that the answer is negative. In the scalar equation (|37|l in 
free space we break (j) into background and perturbation 4> — + Scj), but ignore perturbations of ga/3 and Hq. For 
convenience we shall call 0b simply cj). To first order in S(p we get [c.f. Eqs. HA2|l - (|A4|l ] 

= {h°'^^ + 2^H°'Hf^)S(p,af3 + --- (156) 
= (/i^''</),^0,,)-i/2/i"'3</,,^ (157) 
^ = d\nfi{y)/d\ny (158) 

where the ellipsis denotes terms with Scj) differentiated only once. We have temporarily assumed that H" is spacelike. 
Using Eq. (|23|) we reexpress p5t)|) in terms of the physical metric: 

[g-2,^~a/3 _ (2 _ e-4'>)ii"it/3 + 2^H°'H'^]S<j),^,3 + • • • = (159) 



1. Quasistatic background 



For a quasistatic background, e.g. a quiescent galaxy, H°' is indeed a purely space vector in coordinates that 
reflect the time symmetry. By (I157|l H"' is normalized to unity w.r.t. metric g^p and to e~^"^ w.r.t. (jap- In a local 
Lorentz frame of Qap at rest w.r.t. to those coordinates and appropriately oriented, a generic iJ" will have the form 
e~'^{0, s, 0, Vl — s^}, with s the cosine of the angle between Hi and the positive x axis. Then in a Lorentz frame 
moving w.r.t. the first one at velocity V in the positive x direction 



i?" = e-^{l - V^)-^'^{-Vs, s, 0, v/(l-s2)(l-y2)} 

In this same frame iX" is given by Eq. (|142|l . 

In the eikonal approximation (c.f. Appendix one replaces in a Lorentz frame Scfi-ap 
derivatives. Again interpreting —Kt as u this gives a generalization of (|143|l . namely 



(160) 



-KaKpScf) and drops first 



= AlU^ + 2{BnKn + Bj_Kj_)lO 



Dk{^ - (1 - V^){k1; 



EkI) 



A 



2e^'^ - (1 + 2^5^)1/2 
1/(26*-^ - 1 - 2is^) 



b = 2T/2p/"^ - (1 + 2^5^) 
E = 1 + 2^(1-5^), 



(161) 

(162) 
(163) 

(164) 
(165) 
(166) 



where Ky is the component in the x direction, k±^ is that in a direction orthogonal to x in the plane spanning the x 
axis and Hi, and Ky is the component orthogonal to that plane (we use vector symbols for components to keep with 
previous notation). 

For V = Q (rest frame of matter) there is nothing to distinguish the x axis from HiS direction, so without restricting 
generality we may set s = 1 and speak jointly of k±^ and Ky as a vector kj^. Then the group speed v = \duj / dK^^/"^ 
turns out to be 



V2 



(1 + 2^)2Kjj + K' 

(1 + 2^)kI + k\ 



1/2 



(167) 



(168) 



From Sec. IIIIEl we compute the logarithmic slope 

e(M) - (/i-l)(M-2)/(3M'-6/i + 4) 

whose graph is shown in Fig.|3| In particular, <^ < ^ in a quasistatic region. In the deep MOND regime /i(y) ~ \fyji 
so ^ « i, while in the high acceleration limit [liy) w 1 so ^ w 0. Consequently, in the deep MOND regime, < e.~'^'^ 
with equality for k\_ = 0. In the Newtonian regime vq = 2~^l'^e~'^'^ for all k. Finally, in the intermediate regime 
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X 



-0 . 1 




FIG. 3: The logarithmic slope ^{^) as relevant for quasistationary systems, < /i < 1, and for cosmology, 

2 < fj, < oo. 



2 i/2g '^'1' < vo < {I + 2^)^/^2 ^/^e with lower and upper equality for = and k,j_ = 0, respectively. 
Summarizing, scalar waves propagate subluminally in the frame in which the matter is at rest, provided, of course, 
(/) > 0. 

Since the vector il" defines a privileged Lorentz frame, the form of the wave equation (|159|l is different in different 
frames. Thus we must check explicitly that the subluminal propagation of scalar waves remains valid in all Lorentz 
frames. Since the analytic expressions for general k are cumbersome, we have done so numerically for small positive 
(j). For small V the group speed starts at the value (|167|l . If < 0, ?; for the + mode rises with increasing V 
approaching unity as V ^ 1. By contrast, H > 0, f at first decreases with increasing V only to reach a minimum 
which can be quite narrow and deep for near unity. Beyond the minimum is a critical V past which the + 

mode with positive /«|| is no longer possible. It is replaced by the — mode with opposite sign of /ty , whose v rises as 

V rises beyond the critical V, approaching imity for V ^ 1. 

In summary, provided > as guaranteed (see Sec. 01 for the vicinity of masses embedded in the cosmolo- 
gies studied in Sec. IVIII no case of superluminal propagation is observed for scalar perturbations on a quasistatic 
background. 

2. Cosmological background 

Consider now propagation of scalar perturbations in FRW cosmology. Here il" remains pointed in the time direction, 
and takes the form (|142() in a local Lorentz frame of the physical metric which moves w.r.t. the matter at velocity 

V in the x direction. Since i?" is now timelike, one must change the sign of the argument of the square root in 
definition H157(l . Definition H158(l then requires a switch in sign of the ^ term in Eq. H156(l . We may evidently write 
(/)_Q, — (with C spacetime dependent). It follows from definition H142|l that H°' = -^211" independent of Using 
all this in the modified wave equation (|159|1 . we obtain in the said Lorentz frame, after an eikonal approximation, a 
dispersion relation of the form H143|) with the coefficients A, B and G modified according to the rule e^"^ (2 + 4^)6^*^. 
Thus in the frame C where the matter is at rest (V^ = 0) we now find the isotropic group speed, c.f. Eq. 1147|) . 

wo = (2 + 40"'/' 6-2"^, (169) 

so that according to Fig.|3| for (/> > 0, -yo never exceeds l/\/2. 

For F > we use the analysis leading to Eqs. (|154|I - H155I) with the substitution e^*^ (2 + 4^6'*'^ to conclude that 
the passage to superluminality is forbidden. Numerical plots disclose a behavior of v{V) very similar to the one for 
tensor waves. For + type modes with /ty < 0, w grows monotonically approaching unity for 1^ — > 1. For /«y > modes 
there is a minimum of v at some high V , the narrower and deeper the larger A mode of this type exists only 

up to a critical V beyond the minimum, and is thereafter taken over by the — type mode whose k\\ is of opposite 
sign, and for which v approaches unity as V 1. 

C. Caveats 

Summing up, propagation of weak perturbations of the tensor, vector or scalar gravitational fields of TeVe S is 
always subluminal with respect to the physical metric. We have checked this in detail only for waves propagating 
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on pure cosmological backgrounds or on quasistatic backgrounds. Furthermore, the analysis looked at perturbations 
of one field while keeping the others "frozen" at their background values. A more advanced analysis would have 
examined propagation of joint tensor-vector-scalar modes. This said, no mechanism is evident for the formation of 
causal loops. This under the condition > which, as we have seen, is widely obeyed in flat-space cosmological 
models and quasistatic systems embedded therein. 
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APPENDIX A: ACAUSALITY IN RELATIVISTIC AQUAL 



This comes about because the wave equation for free propagation of ip deriving from the £^ in Eq. ^ (covariant 
derivatives are w.r.t. gafi), 



[/'(LV>„.V',.)ff"''t/3];« = 0, 
leads to the following linear equation for propagation of small perturbations Sip on the background {5"'^, ^Pb}- 

e = dln/'(2/)/dlny 



(Al) 

(A2) 
(A3) 
(A4) 



In Eq. ljA2|l the ellipsis stands for terms where 5^ is differentiated only once. 

For a static background X°' is a unit purely space vector X. In an appropriately oriented Cartesian coordinate 
system in a local Lorentz frame, it will point in the x direction. In such frame Eq. (|A2|) takes the form 



= -5i}M + (1 + '2i)H\xX + S^J^yy + 5lp,,, + 



(A5) 



In the eikonal approximation appropriate for short wavelengths, one sets ip — Ae'*^ and neglects terms with derivatives 
of A or of fcct = ((5.Q. Then Eq. ljA5|l gives 



LO = -kt = [(1 + 2S)k^'^ + ky^ + fc/]^/^ 
The group speed Wg — |9cj/9k|^/^ turns out to be 



(l + 20fcx^ 



ky"^ 



k. 



1/2 



(A6) 



(A7) 



In the deep MOND regime [f{y) = §2/ ' ], 2^ = 1 while in the high acceleration limit [/(y) ~ y], ^ ~ 0. Thus 
whatever the choice of /, < ^ < 1 over some range of y (acceleration). There Ug > 1 if k is not exactly orthogonal to 
X (distances and times measured w.r.t. metric gap)- On the other hand, light waves travel on light cones of gap while 
metric waves do so on null cones of gap- The two metrics are conformally related so their null cones coincide: light 
and metric waves travel with unit speed. Thus most ip waves are superluminal, in violation of the causality principle 
[see SecHlBl. 



APPENDIX B: PROBLEMS FOR PCG IN SOLAR SYSTEM 



The permissible ranges of 77 and e are strongly constrained by the solar system. It can be shown 01 that the 1/r 
force in Eq. H17fl causes the Kepler "constant" of planetary orbits with periods P and semimajor axes a to vary slightly 
with a: 

47r2aVp2 ^ gMq{1 + aoh/Kri). (Bl) 
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Assuming Mq <C Mc, we get k = j so that as we pass from planet to planet, the "constant" varies by a fraction 
~ 2 X /rj. The inner planet periods P are known to better than one part in 10*. Thus rj > 2 x 10^^. 

A stronger constraint comes the perihelia precessions of the planets. The anomalous force in Eq. (|17() generates 
an extra precession Q which in Mercury's case (excentricity 0.206 and a = 6 x IQ-'^^ cm) amounts to 3 x 10~*?7~^ 
arcsec/century. With ij = 2 x 10~^ this already amounts to 0.35% of the Einstein precession, which is measured to 
about that accuracy. Trying to assume Mq > Mc just aggravates the problem. And we are not at liberty to raise rj 
further because for fixed ag, Mc scales as rj . Thus, for example, with r; = 3 x lO""^, the MOND hmit of PCG would 
not apply to galaxies with M < 8 x 10^, a range including many dwarf spirals with missing mass problems ! Hence 
the perihelion precession marginally rules out PCG with a sextic potential. 



APPENDIX C: RELATION BETWEEN DETERMINANTS g AND g 
From Eqs. I|22I23|1 it follows that 

r^gua = e^* + (1 - e-4*)il^il„] (CI) 

Viewing this as multiplication of two matrices, we take the determinant: 

g-^g^e^'*'DetlC{(j),U); IC{(p,U) = I + {1 - e-^'f'p (C2) 

where X is the unit matrix whose components are S'^a while U is a matrix with components il^ilo,. Now both g and 
g are scalar densities, so that their ratio must be a true scalar. Hence Det IC{4>,U) is a scalar. 

In a local Lorentz frame in which the unit timelike vector U" has components {1,0,0,0}, U^s only nonvanishing 
component is U^o = -1. Therefore, Det/C = [1 - (1 - e.-'^'t')] x 1 x 1 x 1 = e"-*"^. Substituting this in Eq. we 
recover Eq. if^ . 



APPENDIX D: RELATIONS BETWEEN m,, nig AND Tg 

To determine Tg one must delve into the region g < R. Assuming that the ideal fluid modeling the matter is at 
rest in the global coordinates, we may write its 4-velocity as Ua = e'^iia — — c'^+'^/^Jq* (see Sec lIIICj) . Let us return 
to Eq. , substitute Tu from Eq. and reorganize the left hand side to obtain 

^-2g.-5?/4(^2^,g,/4y ^ _8^G<P (Dl) 

«P = pe''(2e"2^ - e^^) + m + e„/87rG (D2) 



Integration gives for g > R 



, /4 _ 2Gmg 



\ r{87TGTu + Qtt)e^'/^-''g^dg (D3) 



R. 



rUg = Att ^e'"^''-''g'dg, (D4) 
Jo 

where the integral in Eq. ljD3|l does not contain p since it extends only outside the fluid. 

How much does the "gravitational mass" mg differ from the scalar mass TOs ? For a star the volume integral of p 
is of order the random kinetic energy, which by the Newtonian virial theorem is of order of the gravitational energy 
~ Girig/R. According to Eqs. (|73|l . 1)741) and H92|l this is also the order of the fractional correction to nis or to uig 
coming from the metric factors and e*^. We have not worked out Ttt or Qu in the interior, but from Eqs (|79|l and 
(|83|l we may estimate that the Tu and Qu/^ttG terms contribute to nig terms of 0{kGms^ / R) and 0{Krg^ /GR), 
respectively. Because we assume small k and K, these last two terms are obviously subdominant contributions. We 
may conclude that nig and TOs differ by a fraction of order Girig/ R which is 10~^ for the solar system. 

Let us now calculate <;'e''^^ at g = R using Eq. ifT^ . and and equate the result to —2Gmg/R^ as stipulated 
by Eq. (|D3)| : 
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For the sun Vg/R ^ Gmg/R ~ 10 ^; we see that rg « 2Gmg with fractional accuracy much better than 10 
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